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Abstract: We study the thermal partition function of quantum field theories on ar- 
bitrary stationary background spacetime, and with arbitrary stationary background 
gauge fields, in the long wavelength expansion. We demonstrate that the equations of 
relativistic hydrodynamics are significantly constrained by the requirement of consis- 
tency with any partition function. In examples at low orders in the derivative expansion 
we demonstrate that these constraints coincide precisely with the equalities between 
hydrodynamical transport coefficients that follow from the local form of the second 
law of thermodynamics. In particular we recover the results of Son and Surowka on 
the chiral magnetic and chiral vorticity flows, starting from a local partition function 
that manifestly reproduces the field theory anomaly, without making any reference to 
an entropy current. We conjecture that the relations between transport coefficients 
that follow from the second law of thermodynamics agree to all orders in the derivative 
expansion with the constraints described in this paper. 
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1. Introduction and Summary 

In this paper we explore the structural constraints imposed on the equations of rela- 
tivistic hydrodynamics by two related physical requirements. First that these equations 
admit a stationary solution on an arbitrarily weakly curved stationary background 
spacetime. Second that the conserved currents (e.g. the stress tensor) on the corre- 
sponding solution follow from an equilibrium partition function. 

Landau-Lifshitz |]I|, and several subsequent authors, have emphasized another 
source of constraints on the equations of hydrodynamics, namely that the equations are 
consistent with a local form of the second law of thermodynamics. As is well known, 
this requirement imposes inequalities on several parameters (like viscosities and con- 
ductivities) that appear in the equations of hydrodynamics. It is perhaps less well 
appreciated that the requirement of local entropy increase also yields equalities relat- 
ing otherwise distinct fluid dynamical parameters, and so reduces the number of free 
parameters that appear in the equations of fluid dynamics (see e.g. |1|, 0, for more re- 
cent work inspired by the AdS/CFT correspondence see e.g. 0, ^, m|, |^, H ||, 0' Hill)- 
In three specific examples we demonstrate below that the equalities obtained from the 
comparison with equilibrium (described in the previous paragraph) agree exactly with 
the equalities between coefficients obtained from the local second law of thermody- 
namics. These results lead us to conjecture that the constraints obtained from these 
two a naively distinct physical requirements infact always coincide. We leave a fuller 
investigation of this conjecture to future work. 



In this paper we nowhere utihze the AdS/CFT correspondence. However our work 
is motivated by the potential utihty of our results in an investigation of the constraints 
imposed by the second law of thermodynamics on higher derivative corrections to 
Einstein's equations |jT2[, via the Fluid - Gravity map of AdS/CFT (||13[, see fl^, |5 



for reviews). 

In the rest of this section we summarize our procedure and results in detail. In 



subsection ^3] below we describe the structure of equilibrium partition functions for 
field theories on stationary spacetimes in an expansion in derivatives of the background 
spacetime metric (and gauge fields). In subsection |L^ we then describe the constraints 
on the equations of relativistic hydrodynamics imposed by the structure of the partition 
functions described in subsection |1 . 1| . In three examples we compare these constraints 
to those obtained from the requirement of entropy increase and find perfect agreement 
in each case. 

1.1 Equilibrium partition functions on weakly curved manifolds 

Consider a relativistically invariant quantum field theory on a manifold with a timelike 
killing vector. By a suitable choice of coordinates, any such manifold may be put in 
the form 

ds^ = -e^'^(^) [dt + ai{x)dx'Y + gij{x)dx'dx^ (1.1) 

where i = 1 . . .p. dt is the killing vector on this manifold, while the coordinates x 
parametrize spatial slices. Here a,ai,gij are smooth functions of coordinates x. 

Let H denote the Hamiltonian that generates translations of the time coordinate 
t. In this subsection we address the following question. What can we say, on general 
symmetry grounds, about the dependence of the the partition function of the system 

Z = Tie'^, (1.2) 

on a, Qij and Oj? In this paper we focus on the long wavelength limit, i.e. on manifolds 
whose curvature length scales are much larger than the 'mean free path' of the thermal 
fluid ^. In this limit the question formulated above may addressed using the techniques 
of effective field theory. In the long wavelength limit the background manifold may be 
thought of as a union of approximately flat patches, in each of which the system is in 
a local fiat space thermal equilibrium at the locally red shifted temperature 

r(x) = e-'To + . . . (1.3) 



^Equilibrium Partition functions special curved manifolds or with particular background gauge 
fields have been studied before in |l^, 17, p^ 



(where Tq is the equihbrium temperature of the system and the . . . represent derivative 
corrections, see below). Consequently the partition function of the system is given by 

In Z = I dPxy^^P(T(x)) + . . . (1.4) 

where P{T) is the thermodynamical function that computes the pressure as a function 
of temperature in flat space. Substituting (|1.3|) into ( |1.4|) we find 

In Z = y d^x^^^PiToe^n + ■■■ (1-5) 

The ... in (|1.5|) denote corrections to In Z in an expansion in derivatives of the 
background metric. At any given order in the derivative expansion these correction 
are determined, by the requirement of diffeomorphism invariance, in terms of a finite 
number of unspecified functions of a. For example, to second order in the derivative 
expansion, p dimensional diffeomorphism invariance and U{1) gauge invariance of the 
Kaluza Klein field a constrain the action to take the form 



logZ = W = -Uj d^x^p^^PiToe-n 

+ J dPx^p{P,{a)R + T^P2{a){dia, - d,a,y + P3{a){Vaf) 



:i.6) 



where -Pi (a), ^2(0") and Psicr) are arbitrary functions. It is possible to demonstrate on 
general grounds that the temperature dependence of these functions is given by 

P,{a) = P,{T,e"') (1.7) 

so that 

\ogZ = W = -Uj d^x^,^^P{Toe-n 

+ f dPx^{PiiToe-'')R + T^P2iToe-nid,a, - d^a^f + P:,{T^e-''){Vaf) 

(1.8) 

The discussion above is easily generalized to the study of a relativistic fiuid which 
possesses a conserved current J^ corresponding to a global U{1) charge. We work 
on the manifold ( p..lD in the presence of a time independent background U{1) gauge 
connection 

A = Ao{x)dx° + Ai{x)dx' (1.9) 

and study the partition function 

Z = Tre ^^ (1.10) 



Later in the paper we present a detailed study of the special case of charged fluid 
dynamics in p = 3 and p = 2 spatial dimensions, at first order in the derivative 
expansion, without imposing the requirement of parity invariance. Let us first consider 
the case p = 3. The requirements of three dimensional diffeomorphism invariance, 
Kaluza Klein gauge invariance, and U{1) gauge invariance upto an anomaly^ (see below) 
force the partition function to take the forni'^ 



* ^ ' * ^ %nv ' anom 



T..1 C'o f ,,, T^Ci f ^ T0C2 f ,, (1-11) 

WL, = -r / AdA + -^—^ / ada + -^ / Ada 



27 2 

^AdA + 
3To 6T0 



W' =-( I ^AdA + ^Ada 



where 74,- 



Aq= Aq + flQ 

Ai = Ai~ ^oOj 



:li2) 



( |1.11[) is written in terms of Ai because A,, unlike Ai^ is Kaluza Klein gauge invariant^. 
W^ in ( |1 ■ 1 1| ) is zero derivative contribution to the partition function, and is the 



patchwise approximation to equilibrium, in the spirit of (|1.5|) . W}^^ is the most general 
diffeomorphism and gauge invariant one derivative correction to W^ . Note that W^ is 
the sum of a Chern Simons term for the connection ^4, a Chern Simons term for the 
connection a and a mixed Chern Simons term in A and a. As usual, gauge invariance 
forces the coefficients Cq, C\ and C^ of these Chern Simons terms to be constants, 
( p. .Ill) is the most general form of the partition function of our system that satisfies 



the requirements of 3 dimensional diffeomorphism invariance and gauge invariance. If 



•^In this paper we only consider the effect of f7(l)'^ anomahes ignoring the effects of . for instance, 
mixed gravity-gauge anomahes. A systematic study of the effect of these anomahes in fluid dynamics 
would require us to extend our analysis of charged fluid dynamics to 2nd order, a task we leave for the 
future (see however section n) . It is possible that Ci above will turn out to be determined in terms of 
such an anomaly coefficient. We thank R. Loganayagam for pointing this out to us. 

■^Our convention is 

i I XdY = / d''x^e'^''X,djYk , \fdY=f d'x^2e''d^Yj . 

^The background data can be taken as gauge field A = (^07-^1) with constant chemical potential 
/xo and temperature Tq. Equivalently we can think of the system to have background gauge field 
B = [Aq + /io,y^i) with no chemical potential. These two are equivalent physical statements as /xo 
can be absorbed in the constant part of ^o- 



we, in addition, impose the requirement of CPT invariance of the underlying four 
dimensional field theory then it turns out that Co = Ci = (see subsection |3.6|) . In 
other words, the requirement of CPT invariance allows only the mixed Chern Simons 
term, setting the 'pure' Chern Simons terms to zero. 

^anom ^^ ^^c part of the effective action that is not gauge invariant under f/(l) 
gauge transformations. ^ Its gauge variation under A^ — ;■ A^ + d^(j){x) is given by 

SW^on. = ^ / d'xV^, * (^ A ^) cPix) (1.13) 

As we explain in much more detail below, this is exactly the variation of the effective 
action predicted by the anomalous conservation equation 

V,> = -^*(^A^) (1.14) 

where J is the gauge invariant ^(1) charge current, and * denotes the Hodge dual. 

Let us now turn to parity violating charged fluid dynamics in p = 2 spatial di- 
mensions. In this case there is no anomaly in the system and the parity odd sector 
is qualitatively much different from its p = 3 spatial dimension counterpart. For this 
system we primarily focus on the parity odd sector upto the first order in derivative 
expansion and the manifestly gauge invariant partition function in this case takes the 
form 

lnZ = >V° + >V, (1.15) 



where 



W = - {a{a, Ao) dA + To /3{ct, Aq) da) . 



;i.i6) 



Where Aq and Ai are defined in ( p..l2|) and a and /3 are arbitrary functions. 



It is straightforward, if tedious, to generalize the form of the partition function 
presented in special examples above to higher orders in the derivative expansion. To 
any given order in the derivative expansion, the dependence of \nZ, on gij, a^, a, Aq 
and Ai is fixed by the requirements of p dimensional diffeomorphism invariance and 
gauge invariance in terms of a finite number of unspecified functions of two variables, 
a and ^40. 

We will now define some terminology that will prove useful in the sequel. Let 
Sg denote the number of independent gauge invariant scalar expressions that one can 



^It is striking that the effect of the anomaly can be captured by a local term in the 3 dimensional 
effective action. Note that W^ cannot be written as the dimensional reduction of a local contribution 
to the 4 dimensional action, in agreement with general expectations. 



construct out of a, a^ (and Aq and Ai in the case that the fluid is charged) at n^^ order 
in the derivative expansion. In a similar manner, v" and t" will denote the number 
of 77,*'* order independent gauge invariant vectors and (traceless symmetric two index) 
tensors formed out of the same quantities. Finally let st^ denote the total number of n*'* 
order scalar expressions that happen to be total derivatives (including the contribution 
of a coefficient function) and so integrate to zero ^ It is clear that at n^^ order in 
the derivative expansion, the equilibrium action In Z depends on s" — st" unknown 
functions of two variables. 

1.2 Constraints on Fluid Dynamics from stationary equilibrium 

1.2.1 Relativistic Hydrodynamics 

In this subsubsection we present a lightening review of the structure of the equations 
of charged relativistic hydrodynamics. The equations of hydrodynamics are simply the 
equations of conservation of the stress tensor and the charge current 

v^r,^ = j-.^>, v^> = -^*(^A^), (1.17) 



where J-" is the field strength of the gauge field A in (|1.9D. These equations constitute 
a closed dynamical system when supplemented with constitutive relations that express 
T^j, and J^ as a function of the fiuid temperature, chemical potential and velocity. 
These constitutive relations are presented in an expansion in derivatives and take the 
form 

T'^- = (e + P)ti^V + P(?^^ + TT^^ r = qu'' + J^.^^, (1.18) 

The pressure P, proper energy density e and proper charge density q are those functions 
of T and /z predicted by fiat space equilibrium thermodynamics, tt^i, refers to the 
sum of all corrections to the stress tensor that are of first or higher order in the 
derivative expansion (the derivatives in question could act either on the T, /i, -u^, or the 
background metric and gauge field g^i, and A^). Similarly J^^^^ refers to corrections to 
the perfect fiuid charge current that depend on atleast one spacetime derivative. Field 
redefinitions of the T n and u'^ may be used to impose p + 2 constraints on n^^ and 
•^diss'^ throughout this paper we will work in the so called Landau Frame in which 

^/%, = 0, ^i^jf^ = (1.19) 

Terms in vr^y and J^^^^ are both graded according to the number of spacetime derivatives 
they contain, i.e. 

vr - 7r(^) + 7r(2) + vr^g^ + • • • 

"^diss '^diss,{l) ^ '^diss,{2) ^ "'diss, (3) ^^ • • • 



^For example, the two derivative scalar Vp/i((7)V^cr is a total derivative for arbitrary h{a) 



where the subscript counts the number of derivatives. 

Symmetry considerations immediately constrain the possible expansions for tt^^ 
and J^^^^ as follows. At any given point in spacetime, the fluid velocity u'^ is a particular 
timelike vector. The value of the velocity breaks the local SO{p, 1) Lorentz symmetry 
of the theory down to the rotational subgroup SO{p). In the Landau frame ( [I.17| ) 7r^,y 
may be decomposed into an SO{p) tensor and SO{p) scalar. J^^^^ is an SO{p) vector. 

In order to parameterize freedom in the equations of hydrodynamics, it is useful to 
define some terminology. Let f}, v^ and Sj respectively denote the number of onshell 
inequivalent tensor, vector and scalar expressions that can be formed out expressions 
made up of a total of n derivatives acting on T, u^, /i, gf^^ and A^. It follows imme- 
diately that the most general symmetry allowed expression for n'f^s is given in terms 
of t"^ + s^ unknown functions of the two variables T and /z. In a similar manner the 
most general expression for the ^^j^^c^), permitted by symmetries, is given in terms of 
v^ unknown functions of the same two variables. 

It turns out that the (t^ + s^ + v^) n^^ order transport coefficients are not all 
independent. The requirement that the hydrodynamical equations are consistent with 
the existence of an entropy current that is of positive divergence in every conceivable 
fluid flow imposes several relationships between these coefficients cutting down the 
number of parameters in these equations; we refer the reader to [|I], ^ |, |T^ , for example, 
for a fuller discussion. We now turn to a description of a simpler physical principal that 
appears predict the same relations between these coefficients. These relations may all 
be constructively determined by comparison of the equations of hydrodynamics with a 
partition function. 

1.2.2 Constraints from stationary equilibrium 

As we have explained in the previous subsubsection, it follows from symmetry consid- 
erations that the equations of charged hydrodynamics, at n^^ order in the derivative 
expansion, are parameterized by f] + v"^ + s^ unknown functions of two variables (or 
t"^ + s^ functions of one variable for uncharged hydrodynamics). We will now argue 
that these functions are not all independent, but instead are determined in terms of a 
smaller number of functions. 

It is easy to verify that the equations of perfect fluid hydrodynamics (hydrody- 
namics at lowest order in the derivative expansion) admit a stationary 'equilibrium' 
solution in the backgrounds ( [1.1| ) and (|L^) given by 

M53)(f) = 6-^^(1, 0,...,0), T(o)(f) = ^oe-^ /i(o)(f) = e-'^Ao (L21) 

As explained above, this is also the equilibrium solution one expects of the fluid on 
intuitive ground. At higher order in the derivative expansion this solution is corrected; 



the corrected solution may be expanded in derivatives 

T = T(o)+T(i) + r(2) + ... (1.22) 

/i = /i(o) + ^(1) + /i(2) + • • • 

where u^^)^ -^(n) ^^"^ /^(^i) ^^^ expressions of n^^ order in derivatives acting on a, Aq, 
aj, Aj and (^jj. What can we say about the form of the corrections M/'^w T(^n) and 
/X(n)? Adopting the notation defined in the last paragraph of the previous subsection, 
symmetries determine the expression for u^. in terms of v'^ as yet unknown functions of 
a and Aq , while T and /x are each determined in terms of s" as yet unknown equations 
of Aq and a. 

The stress tensor and charge current in equilibrium are given by plugging ( |1.22|) 



into (|1.2CI|) . The result is an expression for n^'^ and J^^^^ written entirely in terms of a, 



Ao, aj, Aj, Qij and their derivatives. 

This expressions for the stress tensor and charge current so obtained depend only 
on a subset of the transport coefficients that appear in the expansion of tt^'' and J^j^^. 
For instance, the expansion of the n*'^ order tensor part of tx^^ has t^ terms in general. 



n 



When evaluated on (|1.21| ), however, this expression reduces to a sum over t" < t^ 
terms. The coefficients of these terms define t" subspace of the fl dimensional set 
of n*"^ order transport coefficients. We refer to this subspace as the subspace of non 
dissipative transport coefficients. 

In this paper we demand that the expressions for the equilibrium stress tensor and 
charge current, obtained as described in the previous paragraph, agree with the corre- 
sponding expressions obtained by differentiating the equilibrium partition function of 



subsection |L1| with respect to the background gauge field and metric. This requirement 
yields a set of t" + 2v^ + Ss" equations^ that completely determine both the n*^ order 
corrections to the equilibrium solutions T„ /i„ and m^ {v^ + 2s" coefficients in all) as 
well as the f^ + v^ + s" non dissipative hydrodynamical transport coefficients. Note 
that the number of variables precisely equals the number of equations. Dissipative 
hydrodynamical transport coefficients are completely unconstrained by this procedure. 
We emphasize that the shifted equilibrium velocities, temperatures and chemical 
potentials obtained from the procedure just described automatically obey the equa- 
tions of hydrodynamics. By construction, the shifted fluid variables, together with 
the constitutive relations determined above yield the stress tensor that follows from 



^The counting goes as follows. The stress tensor decomposes into one SO{p), tensor, one vector and 
two scalars. The charge current decomposes into a vector and a scalar. Equating the hydrodynamical 
equilibrium stress tensor and charge current to the expressions obtained by varying the equilibrium 
yields 3s" + 2v" + t" equations. 



the functional variation of an equilibrium partition function, and the stress tensor ob- 
tained from the variation of any diffeomorphically invariant functional is automatically 
conserved. Very similar remarks apply to the charge current. 

Let us summarize. In general n^^ and J^^^,^ are expanded in terms of t" + s" 
and v^ transport coefficients, each of which is a function of temperature and chemical 
potential. However fl — f^ + s^r — s" of these coefficients in -k^^ and v^ — w" of these 
coefficients in J^^^^ evaluate to zero on the 'equilibrium' configuration ( p..21|) . The 



remaining t'^ + v'^ + s'^ non dissipative transport coefficients multiply expressions that do 
not vanish on (|1.21| ). Comparison with the equilibrium partition function algebraically 



determines all non dissipative transport coefficients in terms of the s" — st"^ functions 
(and derivatives thereof) that appear as coefficients in the derivative expansion of the 
partition function. In other words the t" + v" + s" non dissipative transport coefficients 
are not all independent; there exist t" + v" + sf^ relations between these coefficients. 

The procedure described above may also be used to derive constraints on the form 
of the fluid entropy current. The entropy current must obey two constraints. First its 
divergence must vanish on all the equilibrium configurations derived above. Second, 
the integral over the entropy density (obtained from the entropy current) must equal 
the thermodynamical entropy that follows from the partition function ( [1.1CI| ). These 



requirements impose constraints on the form of the (non dissipative) part of the most 
general symmetry allowed hydrodynamical entropy current. 

We have implemented the procedure described above in detail in three separate 
examples which we describe in more detail immediately below. In each case we have 
obtained detailed expressions for all non dissipative hydrodynamical coefficients in 
terms of the parameters that appear in the action. In each case, the relations obtained 
between non dissipative transport coefficients, after eliminating the action parameters, 
agree exactly with the relations obtained between the same quantities by previous 
investigations based on the study study of the second law of thermodynamics. 

In the case of parity violating first order fluid dynamics in 3 + 1 dimensions, the 
results for transport coefficients computed from ( p..ll|) match perfectly^ with those of 
Son and Surowka ^ (generalized in [^,0) once we impose the additional requirement 
of CPT invariance. ^ 



^See [^9[ ^ for an alternate (using quantum kinetic theory) derivation of hydrodynamic coefEcients 
related to the chiral anomaly without making any reference to an entropy current. 

^Before imposing the requirement of CPT invariance, we have an additional one parameter freedom 
that is not captured by the the generalized Son-Surowka analysis. The reason for this is that Son and 
Surowka (and subsequent authors) assumed that the entropy current was necessarily gauge invariant. 
This does not seem to us to be physically necessary. It seems to us that an entropy current whose 
divergence is gauge invariant - and whose integral over a compact manifold in equilibrium is gauge 
invariant - is perfectly acceptable. As we explain below, it is easy to find a one parameter generalization 



In the case of parity preserving fluid dynamics in 3+1 dimensions, the results 

obtained from the partition function (|1.6| ) agree perfectly with those of Bhattacharyya 

ID| . Finally, in the case of parity non preserving charged fluid dynamics in 2+1 



dimensions, the results from section ^ agree perfectly with those of |jT8| 



In ending this introduction let us note the following. As we have described at 
the beginning of the introduction, the physical principles that yield constraints on the 
transport relations of fluid dynamics are twofold. First, that these equations are consis- 
tent with the existence of a stationary solution in every background of the form ( |1 . 1|) , 
(11.91). Second, that the stress tensor and charge current evaluated on this equilibrium 
configuration obeys the integrability constraints that follow if these expressions can be 
obtained by differentiating a partition function. In the presentation described above we 
have mixed these two conditions together (as the partition function is the starting point 
of our discussion). However it is also possible to separate these two conditions. For 
each of the three examples discussed above, in Appendix ^ and |C| we present a de- 
tailed study of the constraints on the equations of fluid dynamics obtained merely from 
the existence of stationary solutions in arbitrary backgrounds of the form ( [1 . 1|) , ( |1.9|) . 
In each case we find that all of the relations between transport coefficients, derived in 
this paper, are implied already by this weaker condition. In these three examples, once 
equilibrium exists, the requirement that it follows from a partition function turns out 
to be automatic. We do not expect this always to be the case. In more complicated 
cases we expect the existence of a partition function to imply further constraints than 
those implied merely by the existence of equilibrium. However we leave the study of 
such effects to future work^°. 

2. Preparatory Material 

In this section we present background material that we will need in the main part 



of the paper. In subsection |2.1| we present some Kaluza Klein reduction formulae for 
metrics of the form (|1 . 1|) . In subsection |2.2| we describe the transformation properties of 
various quantities of interest under Kaluza Klein gauge transformations. In subsection 



2.3 we discuss how the stress tensor and charge current of our system is related to the 



of the Son-Surowka solution that meets these conditions, and that gives rise to the additional term Co 



in the partition function (1.11). However it turns out that the requirement of CPT invariance sets Co 



(along with Ci) to zero in ( 1.11 ), so this possible ambiguity is never realized in the hydrodynamical 
description of a quantum field theory. 

"'^'' After this paper was completed we were informed of two upcoming works 

(1) "Towards hydrodynamics without an entropy current" of K. Jensen, M. Kaminski, P. Kovtun, A. 
Ritz, R. Meyer and A. Yarom 

(2) "Triangle Anomalies, Thermodynamics, and Hydrodynamics" of K. Jensen 
which have some overlap with our current work. 



partition function. We also discuss the thermodynamical energy, entropy and entropy 
of our system, and compare these quantities to those obtained from integrals over local 
currents. In subsection P]l| we discuss the relation between consistent currents (those 
obtained from the variation of an action) and gauge invariant currents in systems with 
a [/(I) anomaly. In p.5| we describe how the equations of perfect fluid hydrodynamics 
may be 'derived' starting from a zero derivative equilibrium partition function. 



2.1 Kaluza Klein Reduction Formulae 

As explained in the introduction, in this paper we study theories on metric and gauge 
fields in the Kaluza Klein form 

A^' = {A%x),A\x)) 
The inverse of this metric is given by 



/ij/ 



-e-2- + a2) 



-a 



where the first row and column refer to time and g"^^ is the inverse of gij. Christoffel 
symbols, F, of the p + 1 dimensional metric are given in terms of those of the p 
dimensional Christoffel symbols F by 



00 



00 



-e^''{a.d)a 



T% = d,a - e"'{a.d)aa, + 



6 Jim^ 



Fjo = e'^g^'{--f,, + d,aa,) 



J2a 



r°. = -a F" + — 



1 e 









km 



ttjdiam + aidja. 



j""m 



+ -g'^^dUe'-a^a,) 

Curvature symbols of the p + 1 dimensional metric (e.g. the Ricci scalar R) are 
in terms of p dimensional curvature data (e.g. the p dimensional Ricci Scalar R) 

^^The definitions we adopt in this paper are 

p cr o po" o per I j^a per pa j^a 



(2.2) 

given 
byii 



(2.3) 



R = R+ -e'^-f - 2(Va)2 - 2VV 
K'^ ^ R,V{uy = WW + V^W + \e'^nf^, 



(2.4) 



where fij = diUj — djtti 

Let us define m^ to be the unit normalized vector in the Killing direction. In 
components 

< = e-'^(l,0,...,0) (2.5) 

Let V/c^" denote the projector orthogonal to -u^ 

^yc^^ = ^^^ + <«K (2.6) 

Explicitly in matrix form 

.^ N /^O 

Let us also define the shear tensor, vorticity and expansion and acceleration of this 
Killing 'velocity' field by 

Qk = V.iij^ = Expansion, a^ = (m/^.V)^^ = Acceleration 



^.. ^ p,ap., I'Yj^Ml+^fMa _ Ql^ ) ^ g^^^^ ^^^^^^ 



(2.7) 



u 



K 



p,.p., ^.M,-^v,MA ^ ^^^^^^^ 



A straightforward computation yields 

< = (2.8) 



e 



2.2 Kaluza Klein gauge transformations 

The form of the metric and gauge fields in ( |2.1D is preserved by p dimensional spatial 
diffeomorphisms together with redefinitions of time of the form 

t' = t + (/,(x), x' = x. (2.9) 



We always use the mostly positive signature. 



Under coordinate changes of the form ( p.9| ) the Kaluza Klein gauge field Oj transforms 
like a connection: 

a- = a,; - did). 



I 



Let us now examine the transformation of p + 1 dimensional tensors under the coor- 
dinate transformations ( |2.9|) . It is not difficult to verify that upper spatial indices and 
lower temporal indices are gauge invariant. So, for instance, if A^j^ is any p + 1 dimen- 
sional two tensor, the p dimensional scalar Aqo, the p dimensional vector Aq and the p 
dimensional tensor A^^ are all Kaluza Klein gauge invariant. On the other hand lower 
spatial indices and upper temporal indices transform under the Kaluza Klein gauge 
transformation (|2.9| ) according to 



vl = v,-d4y^, {v'f = v'' + d,(t)V\ (2.10) 

Note that the p dimensional oneforms 

gijV^ = Vi- aiVo 
are gauge invariant. In the sequel we will make heave use of the p dimensional oneforms 

A, = A-a,AQ (2.11) 

This oneform is Kaluza Klein gauge invariant and transform as connections under f/(l) 
gauge transformations. This is the reason that the partition function ( |1.11| ) was written 
in terms of Ai rather than Ai. 

2.3 Stress Tensor and ^(1) current 

The p + 1 dimensional tensors that will be of most interest to us in this paper are 
the stress tensor, the charge current and the entropy current. The stress tensor and 
charge current are defined in terms of variation of the action with respect to the higher 
dimensional metric and gauge field according to the formulas 

5S = j dx^^'y^^^i fAr^jgf^-^ + jf^sA^ (2.12) 

As we have described in the introduction, in this paper we will be interested in the 
partition function In Z of our system on the background ( ^.ip . This partition function 
may be thought of as the Euclidean action of our system on the metric (|1 . 1|) with 
coordinate time t compactified on a circle of length ^. The change of InZ under time 
independent variations of the metric and gauge field is thus given by 



:lnZ = /"rfx^+iy^^ fAr^Jg^" + J^6A,}j 



(2.13) 



It follows that 

(51nZ 

-L nu — —2lt 






The formulae ( |2.14|) are not written in the most useful form for the purposes of this 
paper. As we have described in the introduction, we find it useful to regard our partition 
function as a functional of 

In Z = l^(e^ Ao, a^ A, 9'\ Tq, /io). (2.15) 

By application of the chain rule to the formulas (|2.13| ) we find 
_ Tpe^- 5W ^, _ To fSW 5W 

^J-3{p+i) (>9 V-9iP+i) oAq V-5'(p+i) oAi 

where, for instance, the derivative w.r.t Aq is taken at constant a, ai, Ai, g^\ Tq and 

2.3.1 Dependence of the partition function on Tq and /io 

From the viewpoint of a Euclidean path integral, the parameter Tq in the partition 
function ( |1.10| ) is the coordinate length of the time circle. Moreover, every quantum 
field of charge q is twisted by the phase g^ as it winds the temporal circle in Euclidean 
space. As usual, such a twist is gauge equivalent to a shift in the ' temporal gauge 
field ^0 ~^ ^0 + fJ-o = Aq ^^ holding Ai fixed. It follows that In Z is a function of 
^0) -4j and /io only in the combination Aq and Ai . The dependence of InZ on Tq 
may be deduced in a similar fashion. The Euclidean time coordinate t' = ITq has unit 
periodicity. When rewritten in terms of t', the metric and gauge field retain the form 
(PD with 

e =7^, fli = fliTo, Ao = — 

It follows from all these considerations that 

W(e^A,a^,A,^*^To,/io) = W(^,^,Toa,,A,^'^■). (2.17) 

-to J-o 

We will never use the function W below; all our explicit formulae will be written in 

terms of the function W. Nonetheless ( p^.l7| ) will allow us to relate thermodynamical 

derivatives w.r.t. To and //o to functional derivatives of the partition function w.r.t. 

background fields. 

^^In this formula Aq is refers to the gauge field in Lorentzian space. Note that /ig is gauge equivalent 
to an imaginary shift of Aq in Euclidean space. 



2.3.2 Conserved charges and entropy 

In this subsubsection we will compute the U{1) charge and energy of our system from 
integrals over the appropriate charge currents, and compare the expressions so obtained 
with thermodynamical formulas. 

The U{1) charge of our system in equilibrium is given by 



Q = j (Fx^/^^iJ^ (2.18) 



where the integral is taken over the p dimensional spatial manifold. Let us now define 
the (conserved) energy of our system. Whenever the divergence of the stress tensor 
vanishes, the current —v^T^ is conserved provided v^ is a killing vector field. We 
cannot directly apply this result to the killing vector field t^'^ = (1, . . . , 0), as the stress 
tensor is our paper is not divergence free in general (see ( p. . 1 7p ) . However it is easily 
verified that the shifted current 

J^ = -T,^-AoJ^ (2.19) 

is conserved in equilibrium. As a consequence we define 

E = J dPx^^^rJ% = j cPx,/^^, (-r° - A J°) (2.20) 



Q and E defined in ( ^.ISp and ( |2.2CI| ) may be shown to be Kaluza Klein gauge invariant. 



For instance, the Kaluza Klein gauge variation of the RHS of ( |2.18|) is given by 

/ (Fx^J-gp+iTdi(t) 
= j (Fx^^^iJ^d^<P (2.21) 

(where we have used the fact that the gauge parameter (j) is independent of t, integrated 
by parts, and used the fact that J^ is a conserved current). The gauge invariance of 
E follows from an almost identical argument. 

We will now demonstrate that the expressions (|2.18|) and ( p. 201 ) agree exactly 
with the thermodynamical definitions of the charge and energy that follow from the 
partition function. In great generality, the charge of any thermodynamical system may 
be obtained from its partition function (|1.10| ) via the thermodynamical formula 






where the partial derivative is taken at constant TQ,AQ,Ai,g^^,ai,a. In the current 
context 

(2.22) 



^dW f ( 5W 



dji. 







'SAJx] 



jdPx,/^^,j' = Q 



where we have used ( |2.17| ), J° = — e~^°"Jo — a.iJ'' (this follows from the fact that 
<^ = gooJ'^ + goi'J^) and explicit expressions for Jq and J* listed in ( p. 161) . Let us 
note that, in the presence of anomaly |1.17| current J'^ is neither gauge invariant nor 
conserved^'^ . 

The thermodynamical energy 

,dW 



T, 



'dTo 



+ fj-oQ 



(where the partial derivative is taken at constant ^o,Ao,Ai,g''^,ai,a). may be pro- 
cessed, in the current context, as 



°9Tn 



To 






6ai 6Aq 
j \r^x [(e-^'^Too + a,7^) - A^f\ 







M,- 



(2.23) 



E 



where we have used ( p.l7|) , the fact that —Tq = e ^"^Too + ttjTJ5 ^^^ ^^^ explicit expres- 



sions for Too and Tq in ( 2.16 )). In summary 



E = Tr 



'on 



+ i^oQ 



Q = To 



dW 



Even in the presence of anomaly one can show that the current Jj^ in |2.19 



(2.24) 



remams 



conserved, where J^ is defined as in p.l4[ Thus, the thermodynamic formula p. 23 
holds for anomalous system as well. 
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One can construct a conserved current which is given by 



We conclude that the conserved charge and energy in our system are given, in 
terms of the partition function, by the usual thermodynamical formulae. It follows 
that the entropy of our system should also be given by the standard statistical formula 

S = -g^ (2.25) 



equating ( p. 251 ) with J dPx^—Qpj^iJ^ 



Later in this paper we obtain constraints on the entropy current of our system by 

's- 

2.4 Consistent and Covariant Anomalies 

^^ In this section we discuss the relationship between the consistent charge current (the 
current obtained by differentiating the partition function w.r.t. the background gauge 
field) and the gauge invariant charge current in arbitrary 3 + 1 dimensional f/(l) gauge 
theories with a U{lY anomaly. Readers who are familiar with the issue of consistent 
and covariant anomalies in quantum field theories can skip this section. The equations 
which will be used later are (p:^,(p:^,(p:^. 

In this paper we will have occasion to study field theories in 4 spacetime dimensions 
whose ^(1) current obeys the anomalous conservation 

24 



V^J^ = -— *(jrAj^) (2.26) 



J'^ in ( |2.26| ) is the so called 'consistent' current defined by J^ = ^-. As all gauge 
fields in this paper are always time independent 

*(jr A J-) = -Se-^'e'^^diAQdjAk (2.27) 

(here e^^^ = -^) so that the anomaly equation may be rewritten as 



V^J^ = ^e-'^e'^'^d^Aod^Ak (2.28) 

^^ It follows that the variation of the action under a gauge transformation is given by 

6S = j v/3^^a^0 = ^ y d^x^^,(t> *{TAT) = -^J d^x^,<Pe'^'d,Aod,Ak 

(2.29) 



We now follow the discussion of Bardeen and Zumino p2| to determine the gauge 
transformation property of J'^. The principle that determines this transformation law 



^''We would like to thank S. Trivedi and S. Wadia for discussfons and on this topic and S. Wadia 
for referring us to |2^ . 

-•^^In order to forestall all possible confusion we list our conventions. J-)xi/ — d^Av — d,yA^ , 
*{F A F) = ef^^^^Tfi^Ta/i where e^^^^ — ,\ . The variation of the gauge field Ap, under a gauge 
transformation is given by bAp, = d^cj). 



is simply that the result of first performing an arbitrary variation of the gauge field 
A^ — )■ 6Afj, and then a gauge transformation generated by Sep must be the same as 
that obtained upon reversing the order of these operations. The variation of the action 
under the first order of operations, to quadratic order in variations, is given by 



fv^JJ^iSA, 



(where SJ^ denotes the variation of the consistent current J'^ under the gauge trans- 
formation 6(j)). The reverse order gives 



C f^JiJ^AJ^) _C_ f^JiJ^AJ^) 



24./ ^ SA^ ^ 2AJ ^ 5A^ 



^A = ^[ v/=^'^^ae"''^'<9;30^75 



Comparing the two expressions it follows that under a gauge transformation 

Sr = ^e'^^^^dp(f)T^s (2.30) 

It follows that the shifted current 

J>' = J^- ^e'^'^^^AJ^^ (2.31) 

6 

is gauge invariant. J^ is the current that is most familiar to most field theorists; for 
instance it is the current whose divergence is computed by the usual triangle diagram 
in standard text books. It follows from ( p.31| ) that the divergence of J'^ is given by 

V^> = -^*(J^AJ^) (2.32) 



Using ( p.27| ), the anomaly equations may be rewritten as 



V„> = Ce"e"''d,AodjA, 

Let us summarize. J^ is the gauge invariant current that we will use in the fluid 
dynamical analysis in our paper. It obeys the anomalous conservation equation ( p.32|) . 



On the other hand the non gauge invariant current J'^ is simply related to the action 
W (it is the functional derivative of W w.r.t. A^). These two currents are related by 

(PI. 

To end this subsection we will now derive the stress tensor conservation equation 
( [1.1 71 ) in the presence of a potential anomalous background gauge field. We start by 



noting that the variation of W under an arbitrary variation of g'^'^ and A^ is given by 

6W= f v/=^ (-^Sg^^'^T,, + J>'6A^ (2.34) 



Let us now choose the variations of the metric and gauge fields to be of the form 
generated by an infinitesimal coordinate transformation, i.e. 

General coordinate invariance (which we assume to be non anomalous) demands that 
5W = in this special case. Plugging these variations into ( p.34|) , setting the LHS to 
zero and integrating by parts yields 

j d^xv/^e^ (v^T^, - .P (V, A - V^A) + V^J^a}\ (2.35) 

Using ( p. 26 ) together with the identity 

Xe^""''-^^-^"/? = -^e^'^'^^A^F^pF^, (2.36) 

we conclude that 

V.TJi = 7,,{,P - ^e'^'^^M^J-^^) = F,,r (2.37) 

Thus the two equation of motion of charged fiuids are given by ( p.26| ),( P37D . 

2.5 Perfect fluid hydrodynamics from the zero derivative partition function 

It is well known (and obvious on physical grounds) that the equations of perfect fiuid 
dynamics are completely determined by the equation of state of the fiuid (i.e, for 
instance, the pressure as a function of temperature and velocity). 

In this section we will 'rederive' the fact that the equations of hydrodynamics, at 
zero derivative order, are determined in terms of a single function of two variables, by 
comparison with the equilibrium partition function on a general background of the form 
( p.. 11) . The the results we obtain in this subsection are obvious on physical grounds. 
However this subsection illustrates the basic idea behind the work out in subsequent 
sections. 

At zero order in the derivative expansion, the most general symmetry allowed 
constitutive relations of fiuid dynamics take the form 

rpi^u _ ^^ ^ V)u^'u'' + Vg^", J" = qu^", (2.38) 

At this stage e, V and q are arbitrary functions of any two thermodynamical fiuid 
variables, e, V and q (which will, of course, eventually turn out to be the fiuid energy 
density, pressure and charge density) are as yet independent and arbitrary functions 
of the temperature and velocity. 

We will now show that e, V and q cannot be independent functions, but are all 
determined in terms of a single 'master' function of two variables. In order to do 



that we note that the most general p dimensional gauge and diffeomorphism invariant 



partition function for our system on (|1.1|) must take the form 



W = \nZ = J d'x^^ ^P {Toe-'', e-^o) (2.39) 

for some function of two variables P (it is convenient to regard P as a function of e""" 
and e~" Aq rather than simply a and ^40 as we will see below). The stress tensor and 
charge current that follows from the partition function ( |2.39| ) are easily evaluated using 
( p. 161) . The results are most simply written once we introduce some notation. Let 



a = e'^To, b = c^'^Aq 

Let Pa denote the partial derivative of P w.r.t its first argument, and Pf, the partial 
derivative of P w.r.t. its second argument. Below, unless otherwise specified, the 
functions P, Pa and Pb will always evaluated at (a, 6), and we will notationally omit 
the dependence of these functions on their arguments. In terms of this notation 

Tfj = 0, f = 0, (2.41) 



Comparing the expression for J* in ( p.38|) with the same quantity in (|2.40|) we conclude 
that 

tiA* = e-'^(l,0,...,0) 

Comparing the other quantities it follows that 

V = P, t = -P + aPa + hPb, q = Pb (2.42) 

In the special case of flat space the variables a and b reduce to the temperature and 
chemical potential. It is clear on physical grounds that V, e and q are functions only 
of local values of thermodynamical variables. Consistency requires us to identify the 
local value of the temperature with a and the local chemical potential with b. The 
function P that appears in the partition is simply the pressure as a function of T and 
/x. Standard thermodynamical identities then allow us to identify e with the energy 
density of the fluid and q with the the charge density of the fluid. 

Let us summarize the net upshot of this analysis. Symmetries determine the form 
of the perfect fluid constitutive relations upto three undetermined functions e, V and q, 
of the temperature and chemical potential. On the other hand the equilibrium partition 
function is given by a single unknown function, P, of two variables. Comparison of 
the partition function with the fluid hydrodynamics allow us to determine V, e and q 
in terms of P; as a bonus we also find expressions for the temperature and chemical 
potential in equilibrium on an arbitrary background of the form (|1.1| ) , ( |1.9|) . 



As the results of this subsection are obvious, and very well known. However a 
similar procedure leads non obvious constraints for higher derivative corrections of the 
fluid constitutive relations, as we now explain. 

3. 3 + 1 dimensional Charged fluid dynamics at flrst order in 
the derivative expansion 

In this section we will derive the constraints imposed on the equations of charged fluid 
dynamics, at first order in the derivative expansion, by comparison with the most 
general equilibrium partition function. 

The final results of this section agree with the slight generalization of Son and 



Surowka 0] presented in [^,|3 as we now explain. 

Recall that ^ argued that the hydrodynamic charge currents in field theories with 
a t/(l)^ anomaly must contain a term proportional to the vorticity and another term 
proportional to the background magnetic field. P] used the principle of entropy increase 
to find a set of differential equations that constrain these coefficients, and determined 
one solution to these differential equations. It was later demonstrated that the most 
general solution to these differential equations is a two parameter generalization of the 



Son Surowka result |^,[§]- The further requirement of CPT invariance disallows one 
of these two additional coefficients. 

As we describe in detail below, our method for determining the hydrodynamical 
expansion starts with the action ( ^.11| ), and then proceeds to determine the coefficients 
terms in the charge current proportional to vorticity and the magnetic field in a purely 
algebraic manner. Nowhere in this procedure do we solve a differential equation, so 
our procedure generates no integration constants. However the starting point of our 
procedure, the partition function (|5.11|) itself, depends on the three constants Cq, 



Ci and C2. As we demonstrate below, Ci and C2 map to the integration constants 
obtained from the differential equations of . The third constant Co is new, and does 
not arise from the analysis of 0. As we explain below, this coefficient corresponds 
to the freedom of adding a U{1) gauge non invariant term to the entropy current, 
subject to the physical requirement that the contribution to entropy production from 
this term is gauge invariant. It turns out, however, that the requirement of CPT 
invariance forces Co to vanish. As a consequence this new term cannot arise in the 
hydrodynamical expansion of any system that obeys the CPT theorem. 

3.1 Equilibrium from Hydrodynamics 

In Table (|1|) we have listed all scalar, vector and tensor expressions that one can 
form out of fluid flelds and background metric and gauge flelds (not necessarily in 



Type 


Data 


Evaluated at equilibrium 
T = Toe--, fi = e-Mo, n^ = < 


Scalars 


V.M 





Vectors 


rp rp u 





Pseudo- Vectors 






Tensors 


' fia' ul3{ 2 3 y / 






Table 1: One derivative fluid data 



Scalars 


None 


Vectors 


d'Ao , av 


Pseudo- Vectors 


e'^'d.Ak , e'^%ak 


Tensors 


None 



Table 2: One derivative background data 



equilibrium) at first order in the derivative expansion. It follows from the listing of 
this table that the most general symmetry allowed one derivative expansion of the 
constitutive relations is given by 



TT' 



fiu 



-QOV^u - r]a^y 



JLs = ^{E^- TV^dau) + aiE^ + asP^'^^T + ^^u^ 



CbB^ 



(3.1) 



where the shear viscosity rj, bulk viscosity (, conductivity a and the remaining possible 
transport coefficients bi, b2, 63 and 64 are arbitrary functions of a and Aq. 

We are interested in the stationary equilibrium solutions of these equations. In 
general, every fluid variable can receive derivative corrections in terms of derivatives of 
the back ground data. The equilibrium temperature, chemical potential and velocity 
of our system to first order is given by, 

T = r(o) + 5T = Tqc"'^ + 6T, /i = fi(Q) + 5fi = e~'^AQ + 5fi, 
u^ = mJ,) + Su^ = e~"(l, 0, 0, 0) + Su^", 

Su'^ is determined in terms of (5-u* (which we would specify in a moment) as follows. 
Since both u^ and -uJ^qs is normalized to (—1), we have 

U(o)^5u^' = ^ 5u^ = -ai5u\ (3.2) 

Thus, the nontrivial part of velocity correction 5u^ is encoded in 5u^. 



Solutions in equilibrium are determined entirely by the background fields a^ Aq^ 
Oj, Ai and g^K In Table (0,||) we have listed all coordinate and gauge invariant one 
derivative scalars, vectors and tensors constructed out of this background data. As 
Table (0,|l]) lists no one derivative scalars, it follows immediately that the equilibrium 
temperature field T[x) = e~'^To and chemical potential field iJ,{x) = e~°"Ao receive no 
corrections at first order in the derivative expansion. The velocity field in equilibrium 
can, however, be corrected. The most general correction to first order is proportional 
to the vectors and pseudo vectors listed in Table (|,|]) and is given by 

5u' = -^e^^V.fc + b2B], + hd'a + hd'Ao (3.3) 

where 



fjk 


= djttk 


- dkttj 


Fjk 


= d,Ak 


- duA, 


A, 


= A- 


ajAo 


B'k 


2 


{F,k + Aof,k 


123 


1 





(3.4) 



The fluid stress tensor evaluated on this equilibrium configuration evaluates to ( 2.40 ) 
corrected by an expression of first order in the derivative expansion. The one derivative 
corrections have two sources. 



The first set of corrections arises from the corrections ( |3.1|) evaluated on the zero 
order equilibrium fluid configuration ( p.. 21 ). ^^ Using Table (|^), we then conclude that 



the change in the stress tensors and charge current due to the modified constitutive 
relations is given by 

6Too = 6T^ = 6 Jo = 6T^ = 

1 1 1 (3-5) 

6X = aie'^d'Ao - asToe-'^aV + -{^bAq - -W^/jk + ^Cse'^'F^k 

The second source of corrections arises from inserting the velocity correction (|3.3|) 
into the zero order (perfect fluid) constitutive relations. At the order at which we work 
these velocity corrections do not modify Tqo, Jq or T*-^. A short calculation shows that 
the modification of the stress tensor and charge corrections due to these corrections 



^^When u^ ex (1,0..., 0) the Landau frame condition employed in this paper sets ttoo = T^oi 
jdiss _ Q Consequently Tqo, Toi and Jg receive no one derivative corrections of this sort. 



takes the form 



6Tr 



00 



6 Jo = 5T^ = 



5T^ = -e-(e + P) [^{hAo - h^e'^)e'^^ f^^ + h^e'^^F.j, - hT^e-'^d'a + h^d'A^ 

r1 / 
5J' 



\ {qb2Ao - lqhe-)e^''f,, + \qhe'^^F^, 



- qhToe-^'d'a + qhd'A, 



(3.6) 



The net change in Tq and J* is given by summing ( p.6| ) and ( |3.5| ) and is given by 



6T^ = -e'^ie + P) [^{hA^ - \he'')e'^^ f,k + \h2e''^F,k - hT^e-'^d'a + h0 A^ 



5T 



\{{iB + qh2)Ao - \{i. + qh)e'^y^'f,, + ^(^ij + qh)e'^'F,^ 



(3.7) 



- {qh + a2)Toe-"aV + (^64 + ai)d'Ao 



3.2 Equilibrium from the Partition Function 

We now turn to the study of the first correction to the perfect fluid equihbrium partition 
function ( p.39|) at first order in the derivative expansion. From the fact that Table (H,|l]) 
lists no gauge invariant scalars, one might be tempted to conclude that the equilibrium 
partition function can have no gauge invariant one derivative corrections. We have 
already explained in the introduction that this is not the case; the three (constant) 
parameter set of Chern Simons terms listed in the third line of ( |1 . 1 1| ) yield perfectly 
local and gauge invariant contributions to the partition function, even though they 
cannot be written as integrals of local gauge invariant expressions. In addition to these 
gauge invariant pieces we need a term in the action that results in its anomalous gauge 



transformation property (|2.29|) . This requirement is precisely met by the term in the 
last line of ([LTT|) . ^^ 



With the action (|1 . 1 1|) in hand it is straightforward to use (|2.16|) to obtain the 



^^In order to see this we first note that the last line of (f.ff ) may be rewritten as 

C f 

The variation of this term under a gauge transformation is given by 



C 



d^x^e'^''cj)d,AodjAk 



(3.8) 



in perfect agreement with (2.29). 



stress tensor and current corresponding to this equilibrium solution. We find 



Too = 0, T^ = 0, 



T^ = e-'^e'^'' 






'-\CAI + 2CoAo + C2) V,A, + (2Ci - ^Al - C2A0) V,a,] 

c c 

-AiWjAk + -A^AiW^ttk 
2 f-Ao + Cq\ VjAk + (-Al + C2 j Vj-afc + -A^VjA 



(3.9) 



Using ( p.31[ ) it follows that 
Jo = 0, 



1 



ji = e-''e'^%CA^ + 2Co)V,^ + {-CAl + C2)V,afc] , 



(3.10) 



3.3 Constraints on Hydrodynamics 

Equating the coefficients of independent terms in the two expressions for Tq (p.7[),( pl9D 
determines the one derivative corrections of the velocity field in equilibrium. We find. 

hi 



^' (^i.3C + Av^Co - AuC2 + 4Ci) , 



e + P'3 
63 = 64 = 0. 



(3.11) 



where p = ^ = ^. 

Equating coefficients of independent terms in J* in equations |3.7| and p.lO| and 
using (|3.11| ) gives 

e^ = Cv^T\l - —^^vT) + T^ [(4//C0 - 2C2) - ^(4z.2Co - AvC^ + 4Ci)] , 
3(e + /^j '- e + -T 

^B = CvT{\ - . \ vT) + T(2Co - ^(2z^Co - C2)), 



2(e + P) 



e + P 



«1 = Ct2 = 



(3.12) 



Let us summarize. We have found that the hydrodynamical charge current and 
stress tensor are given by 






(3.13) 



In ( |3.13| ) the viscosities C and rj together with the conductivity a are all dissipative 
parameters. These parameters multiply expressions that vanish in equilibrium and 



are completely unconstrained by the analysis of this subsection. On the other hand 
Cuj and Cb - together with ai and 02 in (lOI ) - are non dissipative parameters. They 
multiply expressions that do not vanish in equilibrium. The analysis of this subsection 



has demonstrated that ai and 0:2 vanish and that (uj and Cb are given by ( 3.13 ). The 
expressions ( p.l3| ) agree exactly with the results of Son and Surowka - based on the 
requirement of positivity of the entropy current - upon setting Cq = Ci = C2 = 0. 



Upon setting Cq = they agree with the generalized results of [^ (see also 



We will return to the role of the additional parameter Cq later in this section. 

3.4 The Entropy Current 

The entropy of our system is given by 

s = ^{n\ogz) 

"^^0 (3.14) 

d^x^e'^'' [CoAiVjAk + SCiT^a^Vjat + 2C2ToA,Vjak] . 

In this subsection we determine the constraints on the entropy current Jg of our 
system from the requirement that ( p.l4|) agree with the local integral 



S= j d^x^^ifg (3-15) 

Notice that the first term in (|3.14|) (the term proportional to Co) cannot be written 
as the integral of a U{1) gauge invariant entropy density. It follows immediately that 
( |3.15| ) and ( p.l4| ) cannot agree unless Jg has a non gauge invariant term proportional 



to Cq. Is it permissible for the entropy current of a system to be non gauge invariant 
(and therefore ambiguous)? Entropy in equilibrium is physical and should be well 
defined. Moreover, if we start a system in equilibrium, kick the system (by turning 
on time dependent background metric and gauge fields) and let it settle back into 
equilibrium, then the difference between the entropy of the initial and final state, is 
also unambiguous. It follows that the entropy production (i.e. divergence of the entropy 
current) as well as ( p.l5| ) are necessarily gauge invariant. However these requirements 
leaves room for the entropy current itself to be gauge dependent. 

Over the next few paragraphs we find it useful to dualize the entropy current to a 
3 form. The addition of an exact form to the entropy three form contributes neither to 
entropy production nor to the total integrated value of the entropy in equilibrium. For 
this reason we regard any two entropy 3-forms that differ by an exact three form as 
equivalent. With this understanding, the unique non gauge invariant entropy 3 form 
whose exterior derivative (the Hodge dual of entropy production) is gauge invariant is 
given by 

AAdA 



The requirement that the exterior derivative of this 3 form to be gauge invariant forces 
its coefficient to be constant. ^^ 

The most general physically allowed form for the entropy current, at one derivative 
order, may then be read off from Table |^ 

J^ = sM^ - uJ^^^^ + DeQu^ + D, (^^ - rp^"9«z/) + DeE>' + Z^.o^ 

+ D^u^ + DbB^' + /le^'^^M.aAX (3.17) 

where /i is a constant 



How is the entropy current ( |3.17| ) constrained by the requirement that its integral agrees 
with (p.l4|) ? The one derivative entropy, as computed from the formula J d^x^/glJg 
has two sources. First, the perfect fluid entropy current su^ has a first derivative piece 
that comes from the one derivative correction of the equilibrium fiuid velocity (see 
above). Second, from the one derivative correction to the entropy current (evaluated 
on the leading order equilibrium fiuid configuration). The terms with coefficients Dq 
and Dc vanish on the leading order equilibrium fiuid configuration. Therefore these 
two coefficients can not be determined by comparing with the total entropy as derived 
from action. All the other correction terms computed from this procedure are parity 
odd, except those multiplying Da and De- It is possible to verify that the integrals of 
the terms multiplying Da and De are nonvanishing and linearly independent. As all 



first derivative entropy corrections in (|3.14|) are parity odd, it follows immediately that 



Da = DE = 0. 

Therefore the zero component of the entropy current at first derivative order is given 
by the following expression. 

JslcorrecUon = sSu' + {-U^B + D b) 5° + {-V^^ + Dj U^ + he^"^'' A^By^A^ (3.18) 



^^Naively, another candidate for a non gauge invariant contribution to the entropy three form is 
given by 

AAd{h{T,ti)U) (3.16) 

where U — it^dx^ and h is an arbitrary function of temperature and chemical potential. But this 
term can be rewritten as follows. 

A Ad {h{T, fi)U)^d {h(T, fi)U AA)- h{T, fi) U A dA 
It follows that this addition is actually equivalent to a gauge invariant addition to the entropy 3 form. 



Using 
To 



u 



—e'^'^aA 



i^j 



ak 



(3.19) 



e^'^'^A^dxA^ = e-^e''^ [Aid^Ak + 2TouaidjAk + Tyaid.ak + di {T^vajA, 



Su^ = —aidu^ = bi 



-y^ttidjak 



- &2 [e'^'aidj {Ak + Touak 



and the expressions for ^b, C^j, ^i and &2 as computed in the previous subsection (see 
( plD and (|3:T^)), we find 



/ 0, X'\/ QA'J s \ correction 



(fxy/gie 



Ajk 



T^ ( 3Ci + /ii^^ + -|^ - z/dfi ) aidjttk 



+ To{2C2 + 2hiy - dB)aidjAk + hAdjAk 



(3.20) 



where 



d, 



T 



Cu^ 



C, 



d.. 



A. 

2^2 



C//3 



2C2U + 2Ci 



(3.21) 



Comparing this expression with ( p.l4| ) we find 

h = Co, dB = 2Coiy, d^ = 2Co//' 
This result agrees precisely with that of Son and Surowka as generalized in |Q 

3.5 Entropy current with non-negative divergence 

In the previous subsection we have determined the entropy current by comparing with 
the total entropy derived from the equilibrium partition function and we have allowed 
for terms which are not gauge invariant provided their divergence is gauge-invariant. 

Now we shall try to constrain the most general entropy current (as given in ( p.l7D ) 
by demanding that its divergence is always non-negative for every possible fiuid fiow, 
consistent with the equations of motion. The analysis will be a small modification of 
because of the new gauge non-invariant term with constant coefficient Co added. 
The steps are as follows. 

• First we have to compute the divergence of the current given in ( p.l7[ ) . The new 
term in the entropy current contributes to the divergence in the following way. 

Cn 



V^ [Coe^''''''A.Vo.Ap] = —e 






~2CoE.B^ 



The full divergence of the entropy current is given by 



+ e{u.V)De + (a.V)D, + (g.V)D, + {E.V)De 
+ De{u.V)Q + Daiy.a) + D.iV.Q) + De{V.E) 



+ 



ODb Db 
dT T 



[B^d^T) + 



da. 



dD 



B 



+ 
+ 



dv 



CTu - 2CnT 



dT 






[u^'d.T) 



dv 



2D 



B 



{uj^d^p) 



2qT 
e + P 

e + P 



uj^^Qn 



D^ + 2DbT 

Db + CTu + 2CoT 



{B,Qn 



(3.22) 



where 



and 



Qf, = df,u ^^ 



E, 
T 



JLs = -^TQ^ + iEE^ + iaa^ + i^uj^ + iBB^, and tt 



iJbU 



-r]a^'' - CQP^" 



As explained in [Q, the divergence computed in (|3.22|) can be non-negative if 
Dg, De, Dc, ^e and ^a are set to zero in the parity even sector. 



Since there is no B^ or u;^ term present in ( 3.22 ), for positivity, in the parity 
odd sector we need all the terms that are lin( 
condition imposes the following 6 constraints. 



odd sector we need all the terms that are linear in 5^ and a;^ to vanish. This 



dD 



B 



D 



B 



dT 
dDB 



0, 



dD,., 2D,. 



dv 



T "' dT 
- CTu - 2CoT = 0, 
2qT 



T 
dD^ 







du 



2Db = 



-e 



B 



e + P 

qT 
e + P' 



D^ + 2DbT = 



Db + CTu + 2CnT = 



(3.23) 



We can determine ^^i, Cb, Db and D,^ by solving these equations. The solution 
is identical to the solution determined from the partition function (as given in 
(|3l2|) and (F3l|)). 



3.6 CPT Invariance 



Field 


C 
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T 


CPT 
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tti 
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— 
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9ij 
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Ao 
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Ai 


— 


— 


— 


— 



Table 3: Action of CPT 



In this subsection we explore the constraints imposed on the partition function ( |1 . 1 1|) by 
the requirement of 4 dimensional CPT invariance. In Table PTB] we list the action of CPT 
on various fields appearing in the partition function. Using this table one can easily 
see that the terms with coefficient Ci and Co change sign under CPT transformation 
while the terms with coefficient C2 and C remains invariant. Thus the requirement of 
CPT invariance of the partition function forces Ci = and Cq = 0. Further it also 
tell us that the function P appearing in the perfect fluid partition function, W^, must 
be an even function of Aq (i.e. that equilibrium does not distinguish between positive 
and negative charges). 



4. Parity odd first order charged fluid dynamics in 2+1 dimen- 



sions 



In this section we will derive the constraints imposed on the equations of 2+1 dimen- 
sional charged fluid dynamics, at first order in the derivative expansion, by comparison 
with the most general equilibrium partition function. The parity even constraints are 
identical to the ones found in 3+1 dimensions (which has been extensively discussed 
in §^. Therefore in this section we shall primarily focus on the parity odd constraints 
which are qualitatively much different from their 3+1 dimensional counterpart. These 
constraints have been obtained using a local form of the second law of thermodynamics 
in 



18|, which we shall reproduce starting from the most general equilibrium partition 



function. 



4.1 Equilibrium from Hydrodynamics 



Partially borrowing some notations from equation (1.2) in |T^, the the most general 



Pseudo-scalars 


e'^diAj , e'Wittj 


Pseudo- vectors 


e'WiAo , e'W'a 


Pseudo-tensors 


None 



Table 4: One derivative parity odd diffeomorphism and gauge invariant background data. 
Here e*-' is defined so that e^^ = —;=. 



symmetry allowed one derivative expansion of the constitutive relations is given by ^^ 

(4.1a) 
(4.1b) 



T^-^ = eu^u" + {P- CVc.m" - XbB - Xn^) P^" - r]a^"' - fja^"' , 



J^' = pu^" + aV^ + aV^ + xeE^ + xtT^ ■ 
The various quantities appearing in the constitutive relations ( |4.1| ) are defined as 






and 






T' 



(4.2a) 
(4.2b) 



(4.2c) 



EP 



^tiup 



UyEp, 



yt^ 



^tiup 



UuVp, 



^M- = - {e'^'^Pu^a/ + e^"^n„a/) , f'^ = e^^^w.VpT. 



(4.2d) 
(4.2e) 



The thermodynamic quantities P, e and p are the values of the pressure, energy density 
and charge density respectively in equilibrium. The transport coefficients xb, Xn, Xe 
and xt are arbitrary functions of a and Aq. The only non-zero quantities in equilibrium 
are B, u, E^ and T^. The rest of the first order quantities appearing on the RHS of 
( [4 .11) vanish on our equilibrium configuration. In Table ^ we list all the parity odd 
diffeomorphism invariant background field data. In Table ^ we list the first order 
quantities occurring in the constitutive relations that are non-zero in equilibrium and 
express them in terms of the background metric and gauge fields^°. 

We are interested in the stationary equilibrium solutions of the fiuid equations aris- 
ing from constitutive relations (|4.1|) . Solutions in equilibrium are determined entirely 
by the background fields a, Aq, ai, Ai and g^^ . Just like in 3+1 dimensions the zeroth 
order solution of the fiuid fields are given by 



u 



{e-^0,0}; T(o) = roe 



/^(o) 



'Ao. 



(4.3) 



^^Note that in this constitutive relation the parity even constraint, namely the Einstein relation, 
have already been taken into account. 
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In the following, we shall use e 



12 



1 

%/S2 



Type 


Data 


Evaluated at equilibrium 


Pseudo-Scalars 


B 


e'^diAj + Aoe'Wittj 




n 


-e'^e^Wittj 


Pseudo- Vectors 




Eo = 0, E' = e-^e'WjAo 
To = 0,f' = -e-^e'^dja 


Pseudo- Tensors 


none 





Table 5: One derivative fluid data which are non-zero in equilibrium. 



The unit normalized vector in the killing direction is, 



w 



K 



:i,o,o) 



In Table(^,^) we have listed all coordinate and gauge invariant one derivative parity 
odd scalars, vectors and tensors constructed out of this background data. Since there 
are 2 pseudo-scalars and 2 pseudo- vectors we can have the following most general parity 
odd corrections to the fluid fields at first order 



u. 



u 



(0) 



T = T< 



(0) 



tb Bx + Tn ^K, 



(4.4) 



fi = /i(o) + mB Bk + mn fix, 



where ^e, i^t, tb, tq, ms and rriQ are taken to be arbitrary functions of a and Aq to be de- 
termined by matching with the equilibrium partitions function in § [4.3| . E'^, T^, Bk, ^k 
are the vectors and scalars, defined in equations f4.2a| and |4.2c| , velocity u replaced by 

Uk- 

Just like in 3+1 dimensions the fluid stress tensor evaluated on this equilibrium 
configuration evaluates to (|2.40|) corrected by an expression of first order in the deriva- 



tive expansion. The one derivative corrections again have two sources. 

The first set of corrections arises from the corrections ( [4.1D evaluated on the zero 
order equilibrium fluid configuration ( |4.3| ). ^^ The second source of corrections arises 
from inserting the fluid field corrections in ( [4.4|) into the zero order (perfect fluid) 
constitutive relations. The net change in the stress tensor and the charge current at 
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When w^ ex (1,0..., 0) the Landau frame condition empioyed in this paper sets ttoo = i^oi 
jdiss _ Q Consequently Tqo, Tq^ and Jg receive no one derivative corrections of this sort. 



first order is obtained by summing these two contributions and is given by 

+ (e + P)Uu^,~^E'^ + u^^o)K) + (e + i^)eT(<o)^^ + «(o)^x)- 



(4.5) 



SJ^ = (^^r^ + ^m^J BkU>1^^ + \^—rn + ^mn j fi^^nf^) 



+ iXE + P^e)E'^ + iXT + P^T)fl 



For future reference it will be convenient to to write down some of the components of 
the stress tensor and current in (|4.5|) purely in terms of the background fields using 
the expressions listed in the third column of Table p. 

,^,, fdP dP \ „ „. fdP dP \^ ,,. 

de de 

ST^ = {-{e + P)CEe'WjAo + (e + P)^Te'W,a) 

Sr = e-- {{XE + piE) e'^d.Ao - {xt + piT)e'^ d,a) . 

(4.6) 



5Too = e^" {{%^B + l^mB^) [e'^d^Aj + A,e'^d,a^) - e"^ T^Tf, + j^m^ ) e^^^.a, ) , 



4.2 Equilibrium from the Partition Function 

We now turn to the study of the first correction to the perfect fluid equilibrium partition 
function ( ^.391 ) at first order in the derivative expansion. From the fact that Table 



lists two gauge invariant Hence the most general parity odd equilibrium partition 
function is given by 

W = i / (a(a, A^)dA + To/3(a, Ao)da) , (4.7) 

where a and /? are two arbitrary functions in terms of which all the 4 transport co- 
efficients and the 6 first order corrections to the velocity, temperature and chemical 
potential are to be determined. 

With the action ([4.7|) in hand it is straightforward to use (|2.16|) to obtain the stress 



tensor and current corresponding to this equilibrium solution. We find 
T^ = 0, 






Too = -Toe'' TT^'^d^Aj + To^e'^d.a, , , 



n - Toe- ( ( To| - Ao|) ^W,a + (to|^ - Ao|^) ^^O^Aj , (^i 



a^o dAa 



Jo = -Toe^ ( TTir^'^d^A^ + To-f e^^a.a, , , 



4.3 Constraints on Hydrodynamics 

In this subsection we shall equate the coefficients of independent terms in ([4.5|) (or 
( [4.6| )) with those in ( [4.8p , to determine the first order transport coefficients and ffuid 
corrections in terms of the two arbitrary functions in the action ([4.7|) . 

The fact that T*-' as evaluated from the action ( [4.7| ) vanishes immediately implies 
from 



dP dP 

dP dP 

Xn = ^m + -^mn, 



(4.9) 



Comparing Too from ( [4.6|) and (|4.8| ) we have 
de dt ^ -„da 



(4.10) 



Comparing Tq from ( [4.6| ) and ( |4.8| ) we have 

Toe- / dp da 



e + P) V dAo M, 

, _ Toe- ( d[5 da\ 

Comparing Jo from ([4.6|) and ([4.8]) we have 



(4.11) 



dp dp da 

-JB + T^m,B = To- 



dT "^ dp " "Mo' 

dp dp _^ ( dp da 

7^T"n + -T^mQ = -Toe To-- Aq- 



dT '' dp '' " V ^^0 ^^0 



(4.12) 



Finally, comparing J* from ( [4.6p and ( ^4.8| ) we have 

^ da 
oAq 

da 
Xt + PKt = -^• 



(4.13) 



In order to compare the constraints obtained in this section with that in [|TB| we 
find the following thermodynamical identities useful 



dP _ fdPdp _dP dp\ I fdp de _dp de\ 
a7~ \dTdp^~dpdf) / \dpdf~dTdp) 



dP 

dp 



dP de dP d€\ / f dp de dp de 
dfdp ^ 'dpdf)/ \dpdf ~ dfdp 



(4.14) 



Now solving for tb,tq, uib and m^ from (|4.1CI|) and ( [4.12|) , plugging the answer in 
to ( |4.9D and using the thermodynamical identities ( [4. 141 ), we have 



a7V 



„da\ dP ( ^ da 

-(7 \ I / rji 



XB = ^|-Toe-'-j+-^T„g^ 



dP frr. 



da da J J op 



Toe-MTol^-^'" 



Mn 



dA, 



(4.15) 



Finally plugging in the values of ^e and ^t from ( ^.11| ) into ( [4. 131 ) we have 

da\ p 



Xe 
Txt 



n 



dAo 

r^ rrda 



P 



Tne" 



dp da 

-'0771 ^0 



Mn 



P 



P 



Toe 



-2a 



dAo 
da\ 



^dp 



(4.16) 



Thus through ( [4.15[ ) and ( |4.16| ) we are able to express the 4 transport coefficients in 
terms two arbitrary functions in the action. The two dimensional manifold of allowed 
transport coefficients is identical to that in equation (1.8) in [TH ^^. In particular it 
easy to eliminate a and /3 from ([4.15|) and ( |4.16| ) so as to obtain the following relation 
between the transport coefficients 



Xb 



e + P 



Xn 



dP dP ^ 

-^-Xe + ^-T XT- 
dp de 



(4.17) 



Note that this relation is identical to equation (4.29) in 118 



^^Note that the additional function fn{T) may be reabsorbed into a redefinition of AIa{ii,T) in 
equation (1.8) in lli. 



4.4 The Entropy Current 

In this system In Z is simply given by the action 



In Z = - / (a(a, Ao)(iA + To^(a, Ao)rfa) (4.18) 



The entropy that follows from this partition function is 

s = 4(r„i,>z) 

1 /■ / da , 5a \ , , _ / 9/3 , 9/3 , , 



(4.19) 



We will now utilize ( [4.19| ) to constrain the hydrodynamical entropy current of the 
system. The entropy current must take the 'canonical' form su^ — vJ^^^g corrected by 
first derivative terms. As in the rest of this section we keep track only of parity odd 
terms. It follows from Table ^ that the most general one derivative entropy current is 
given by ^^ 

Jf,) = su^-^ [xeE^ + XtT^^ + [tieE^ + nrf^^ + {ubB + n^fi) u^ . (4.21) 

nE^TiT^TLB and n^ are functions of temperature and the chemical potential. On sub- 
stituting the equilibrium values of temperature and chemical potential they turn into 
functions of a and A^. We find it convenient to define the quantities 



nE = nE- —Xe + s^E, 
nT = nT- —XT + s4r- 
in terms of which the first order part of the entropy current is given by 



(4.22) 



(/ Qg ()g \ t Qg Qg \ \ 

+ ueW + firf^ 

As we have explained above, the entropy current is necessarily divergence free in 
equilibrium. This condition yields one condition 



^■^The map between the corrections to the entropy current in our paper to that in p8| , considering 
first order terms which are non-zero in the equihbrium, is given by 

m = v\ + — ; n-B = i>2 + 1^4 + -^\ nB = Vi] riQ = v^. (4.20) 



( [4. 241 ) is solved by the ansatz 



where tt, is a arbitrary function of a and ^40. Plugging in this solution, we now have 
a 3 parameter set of entropy currents parameterized by ub, n^j and n. The entropy 
( [4. 191) is an integral over the two parity odd scalars of the system. Equating ( |4.19| ) with 



J (PxyJ—g^Jg, and equating the coefficients of these two scalars, yields two equations 
for Ub, Ui^ and n. We now explain how this works in more detail 
Using the fact 

E° = -e-'^e'^ttidjAo; f° = e'^'e'^aidja, (4.26) 

and the expressions of B and Q in terms of the background field (from Table ^), the 
entropy can be evaluated from the entropy current in a manifestly Kaluza-Klein gauge 
invariant way 



S = j Sx^^.Jl 



Os ds \ 

■^"Tb + -K-f^B + riB j {dA + Aoda) ^^^ 27) 

ds Os 

-^m + o-"^f^ + '^^ 1 ^''^^ - T^nda | . 




Comparing ( |4.27| ) with ( [4. 191 ) and using the thermodynamic identities 

ds \ ds u , , 

we get the following simple expressions 

Ub = a 

(4.29) 
Tin = Toe-'' [Aoa - 2/3 - n) 

In other words, we have managed to evaluate ub, and one linear combination of n^ 
and n in terms of the functions, a and /3, that appear in the partition function of our 
system. Note that we have not been able to completely determine the non dissipative 
part of the entropy current using our method (the method based on positivity of the 
entropy current achieves this determination). However, it straightforward to verify 
that the constraints (equations 3.11, 3.17, 3.18, and 3.20) in [|^ on the corrections to 
the entropy current from the second law of thermodynamics, are consistent with the 
relations (^^) and ( ^125]) . 



4.5 Comparison with Jensen et.al. 

In this subsection we shall give a precise connection between partition function coeffi 
cients a, (3 in equation 



and Mq, M that appears in [|T8|. Comparing equations 
( [4.15D ,( ^?TBD with equation 1.8 of []T^, we get the following differential equations 



dM 



Tn 



da 



BM dM 



dT 



dfj, 

dMn 






-M 



-T I r,|f - A, ^" 



T 



OMn 
dT 



on 



Te-" Tr 



dAo 
dp 



'da 



A< 



dA, 
da 



(4.30) 



By solving first two equations in [4.30| , we get 

M 



a 



+ c 



(4.31) 



Toe- 
where c is some constant. Infact the entropy current presented in this paper matches 
that of |]TB| only if we set c = (see equation (3.22) of []T3). Solving last two equations 



in [4.301 , we get 

where ci is some other constant. 



Tf^Mn + -, / UT,e-)e''^da + ^ ——^ + ci, (4.32) 



Also comparing ( |4.15| ) and ([4.16|) with equations (3.17) and (3.18) in [^ one can 
express the entropy current corrections in terms of a and (5 in the following way 



T- 



rp2 



dv5 
dT 

dh 



,di)i 
di>4 



+ Vl 



-Toe-'' 



da 
d^' 



Tn 



da 



dAn' 



Toe 



-2ct 



Tc 



d/3 
d^ 



A, 



da 
d< 



'a 



(4.33) 



^ „ , ^ dB , da 
+ z>3 = -Toe-'' To-^ - Ao- 



di§) " ^ V dAo ""dAo, 

Note that with this identification, the equation (3.20) in [^ automatically follows. 

4.6 Constraints from CPT invariance 

Imposing CPT invariance of the partition function [4.7| constrains the form of the other- 
wise arbitrary functions a, j3. (and hence all transport coefficients determined in terms 
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Table 6: Action of CPT 

of a and /3). Note that we define parity in 2+1 dimensions as Xi — )■ —Xi and X2 — )■ X2- 
In Table ^ we list the action of CPT on various fields appearing in the partition func- 
tion [4.7| . Based on the Table ^ we see, in ^^ "rfA" changes sign where as "da" does 
not, which implies a is odd under CPT and /3 is even under CPT. 



5. 3 + 1 dimensional uncharged fluid dynamics at second order 
in the derivative expansion 

In this section we will derive the constraints imposed on the equations of uncharged 
fluid dynamics, at second order in the derivative expansion, by comparison with the 
most general equilibrium partition function. We do not assume that our system enjoys 
invariance under parity transformations. We do, however, assume that the fluid enjoys 
invariance under CPT transformations. 

Before getting into the details let us summarize our results. Symmetry consid- 
erations determine the expansion of the hydrodynamical stress tensor upto 15 parity 
even and 5 parity odd transport coefficients. It turns out the 7 of the parity even and 
2 of the odd terms vanish in equilibrium. In other words, on symmetry grounds our 
system has 7 parity odd and 2 parity even dissipative coefficients. In addition we have 
8 parity even and 3 parity odd non dissipative coefficients. The most general second 
order fluid dynamical partition function, on the other hand, is given in terms of three 
functions of a ^^. It turns out that this partition function is automatically even under 
parity transformations. As a consequence, implementing the procedure spelt out in the 
introduction, we are able to show that the three nondissipative parity odd coefficients 
all vanish. In addition the 8 nondissipative parity even coefficients are all determined 
in term of three functions. In other words we are able to derive 5 relations between 
these 8 parity even coefficients. 
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The only possible first order contribution to the partition function is the term proportional to Ci 



in (1.11). As explained in subsection 3.6 the requirement of CPT invariance forces Ci to vanish. It 
follows that there are no first order contributions to the partition function for an uncharged system. 
We thank S. Dutta for discussions on this topic. 



The problem of constraining fluid dynamics at second order in the derivative ex- 
pansion, using the principle of entropy increase, was studied by one of the authors of 
this paper in jTO]. In that work the fluid was assumed to enjoy invariance under parity 
transformations. It was demonstrated that the principle of entropy increase indeed 
implies 5 relations between the 8 non dissipative transport coefficients. It turns out 
that the five relations determined in this paper agree exactly with those of [p!0| . 

Even from a practical point of view the method used in this paper appears to have 
some advantages over the more traditional entropy method utilized in |T0[. To start 
with the algebra required for the analysis in this paper is considerably less formidable 
than that employed in |10|. As a consequence we are able, rather effortlessly, to gener- 
alize our results to allow for the possibility of parity violation. Such a generalization 



would involve considerable extra effort using the method of [jTO[, and has not yet been 
done. 



5.1 Equilibrium from Hydrodynamics 



In Tables 1, 2, 3, 7 of |T0[, all scalar, vector and tensor expressions that one can form 
out of fluid flelds and background metric (not necessarily in equilibrium) at second 
order in the derivative expansion are listed. It follows from the listing of these tables 
that the most general symmetry allowed two derivative expansion of the constitutive 
relations is given by 



n 



flU 



- v(^f^u - CPf^uQ 



+ T 



T (u.V)cr(^.,.) + HiR{f,u) + K2K{^u) + Aq Q(T, 



/ii/ 



+ Al CT/ "-aaiy) + A2 CT/ "-iOaiy) + A3 U/^^-Uau) + A4 tt/^a 



'{fi ^au) 



{fi ^av) 



^{tx'^v) 



+ TP, 



+ T 



fii/ 
4 



(5.1) 



Ci(M.V)e + C2R + Cs^oo + 60' + 6^' + C-sco^ + C^a" 



^5,t»+<55^M<^a„r 



>- i=l 



Type 


Data 


Evaluated on equilibrium 


Pseudo-Scalars 


l^ap 


le'^e^^'d.afjk 


Pseudo- Vectors 


u.Vlp 






Pseudo- Tensors 


f^^^ - 1 a 
^(2) ^ e^P'^^uxapao^^^g^y^, 
tfu = e'-P^^uxVpOo^^^g^yp, 

ttMU = UbR^pei.ycdqU'^ 







Table 7: Two derivative parity violating fluid data(Here di^2 are function of a determined 
by evaluating t^^, on equilibrium, but we will not need there explicit expression.) 



where 



uP = The normalised four velocity of the fluid 
pp'^ = gp'^ + uPu'^ = Projector perpendicular to uP 
= V.u = Expansion, a^ = (u.V)m^ = Acceleration 

^M. ^ p,^apu„ , ■— ^ ' ^/3W^ _ ® g^^^ ) = Shear tensor 



2 3^"'^ 



(5.2) 



^^lu _ pfiapp 



^0 (VaUp - VpUc 



Vorticity 



RP'' = Rp'"'\aUb, Rp" = R^P^'^gab (R"^""^ = Riemann tensor) 



a' = ap,aP'', u^ = u^^u'^p 



and 



A 



{nu) 



papp 



Aa(3 + A 



13a 



A„.P 



ab 



gap For any tensor A 



p.U 



The parity odd terms in the last bracket in ( ^.1| ) are deflned in Table |^. 

The expansion (|5.1|) is given in terms of 15 undetermined parity even and flve 
undetermined parity odd transport coefficients, each of which is, as yet, an arbitrary 
function of temperature). 

We are interested in the stationary equilibrium solutions of these equations. Solu- 
tions in equilibrium are determined entirely by the background fields a^ Ui and g^^ . In 
Table (^,|ID we have seen that the 6 and a pi, evaluates to zero in equilibrium. This sets 



seven of the fifteen parity even terms in equation |5.1| to zero. Two of the five parity 
odd terms two terms {tp^ in table 0) evaluate to zero in equilibrium. The remaining 
8 parity even and 3 parity odd coefficients are non dissipative; the non dissipative part 



of n^i, is given by 

n 



/ll^ 



l^lR{nu) + l^2K{^y) + X^UJi^^'^Uau) + A40(^ay) 



T 

+ PAC2R + c^R^o{uy + 6^^' + ^40') 

+ 5iti, + 54% + ^s^M-a^,/" (5.3) 

In order to proceed further, we list all coordinate invariant two derivative scalars, 
vectors and tensors constructed out of background data are listed in table (H). The 
temperature and velocity in equilibrium receives correction at second order. The most 
general symmetry allowed form of corrected temperature and velocity is 

\m=l / 

\m=l / 

where, VmiV) and S'j(S')are Vectors (pseudo) and scalars (pseudo) respectively that are 
listed in table |^. Also 60 can be fixed following equation ^]2| as, 

60 = 1 - e^a. Y^ v^Vi^^) + vV J (5.4) 

\m=l / 

As in previous sections, the stress tensor in equilibrium received corrections at 
second order in the derivative expansion. The two derivative corrections have two 
sources. The first set of corrections arises from the corrections ( ^.1|) evaluated on the 
zero order equilibrium fluid configuration. Using 

u'= = -^^r, a^ = g'^dm(y, (5.5) 



R<ij> = Rij - ViaWja - WiVja + -/^ ^ fjke'^" 
-^(i?-(Va)^-VV+i/V-),., 

a<iaj> = VicrVjCr - -{Wafgij (5.6) 



Scalars 


S, = R, S2 = V'a, Ss = {Vaf, S, = pe^^ 


Pseudo- Scalar s 


S = e'^%afjk 


Vectors 


V, = e-V^aP, V2 = e<^V,r , 


Pseudo- Vectors 


V^ = U.hie^^' 


Tensors 


RijJi''fkj,ViVj(7,V^aVja 


P seudo- Tensors 


di^iaej)kif^\ V{iej)kif^ 



Table 8: Two derivative background data 



we find that these corrections are given by 



n 



eg 



ai i Rij - —Qij ) + ^2 ( ViVjCr - V^agij j + 03 ( ViaVjCX - 
+ «4 (ftfkj + ^g^^ e"' + gi, Ur + &2 VV + bsiVaf + bj 



(Va) 



-9ij 



(5.7) 



+ ^T{6, + 6di)e''d^aej)kif^' + \Td25,e''tki^^^^f^' + hrS.e'' g^.e^ikd'^af^ where, 



^ = C2+gKl, ^ = 2(/«2-Ki) -2C2 + C3 
-^ = 77(«;2 - Ki + A4) - 2C2 + CS + ^4 

1 o 



T 



l(A3-K2-2^0 + ^(C2 + C3-e3), Y 



Ki 



— = K2-Ki, — = K2-A;i + A4, — = -(A3-2/ti -k;2). 



(5.8) 



Here, the indicies are contracted with the lower dimensional metric gij and its inverse. 
The coefficients are determined by evaluating the tliu in equilibrium, but we will not 
need the explicit expressions. 

The second source of corrections arises from inserting the velocity correction ( p.3| ) 
into the zero order (perfect fluid) constitutive relations. We find that the modification 
of the stress tensor due to these corrections is given by 



(5.9) 



Too = Tq ( y ^ traSjn + to 

The net change in Tq and J* is given by summing ( p. 7|) and ( |5.9| ) and is given by 

T' = Prg'' (J2 ^™^- + ^^) + K 

Too = T^^ (j2^mSm + is) (5.10) 



where Wj^ was listed in (|5.7|). 

5.2 Equilibrium from the Partition Function 

We now turn to the study of the first correction to the perfect fluid equilibrium partition 
function ( p.39| ) at second order in the derivative expansion. We observe that the Table 
(^ lists four scalars and one pseudo-scalar. The most generic partition function for 
this system at two derivative order is, 

W = logZ = -^J d\ v/^ [A(Toe-'^)i? + T^p2{Toe-^)f\,r + P^{T,e-''){da) 

where PiiTQe'") = Pi{a) and P'^ = — ^ (2 = 1,2,3) 

da 

(5.11) 

where Pi, P2, P3 are three arbitrary function of a and from now on we will remove 
the explicit dependence. In partition function, the fourth scalar V^a and the pseudo- 
scalar eijkd^af^^ do not appear as they are total derivatives. 

With the action (|5.11 ) in hand it is straightforward to use analog of ( 2.16|) for 



uncharged case to obtain the equilibrium stress tensor. We find 

T^ = TPi {R^ - ]^Rg'^) + 2T^TP2 {f'f.k - \fg'') + T^P^ - P^) (VV VV 
- \iy<yf9'') - TP[{V'V^a - g'^V'a) + h'P'liyafg'^ 

Too = ^{PiR + T^P^f - P-^iVaf -2PsV'a)) 

T^ = 2T^T{P^V,aP' + P^V.r) , (5.13) 

where ' denotes derivative with respect to a. 

5.3 Constraints on Hydrodynamics 



Comparing non trivial components of the stress tensor Tq, Tqo in equations |5.1 0| , |5yT3| and 
equating coefficients of independent sources one obtains the velocity and temperature 
corrections in terms of the coefficients P appearing in |5.11| . We find 



Vl = ^-^2> ^2 = ^--^2, V = 0, 

1 , 1 1 , T^ , ~ 

^1 = 77^— Ai h = --^—Ps, h = - F3, ^4 = 77^—^25 ^ = 0. 



(5.14) 



^^The stress tensor can be evaluated as 



^ Tpe^- SW ^, _ To SW 

-too — / -; — , -^o ^ , =-; — I 

V'~fi'(p+i) "'^ V^-9(p+i) "°* 

T- = -^^=9''9^-^. (5.12) 



Now comparing Tij in equations p.l0| , [5.13| , and using expressions for temperature 
corrections, one can express the transport coefficients in terms of the three coefficients 
P appearing in |5.11| . We find 

TPi, a2 = -TPi a^ = -2T^P2, a^ = T{Ps - P['), 

P3, 64 = -777^^2, (5.15) 
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h 



2P 



TT 



-P[. b2 



-TP[' 



P. 



T 



P' 



Ptt 
Si 



Sd = Sk 



= 0. 



2 2Pj"j' 

One can ehminate the coefficients P' s from above set of relations which gives five 
relations among transport coefficients, 



ai + 02 
1 Ptt 



TdTtti = 0, 



TPtt , 1 / 



Pt 



h + (02 - TdTa2) - as = 0, 4 



TdTtti) 
I Ptt 







Pt 



bd — (3a4 — TdTCLd) 



' ' hi) (02 - r5r«2) - TdT{a2 - ^(9^0.2) - (0.3 



-&3 + (- 



0, 
■ T^tOs) 



= 0. 

(5.16) 



Note that parity odd contributions, both to the equilibrium value of the temper- 
ature and velocity, as well as to the constitutive relations, are forced to vanish. The 
simple reason for this is that the most general two derivative correction to the partition 



function ^.ll] is parity even. Note also that the eight parity even non dissipative trans- 
port coefficients are all determined in terms of the three functions that parameterize 
the two derivative partition function. This leaves us five relations among the transport 
coefficients; these relations may be obtained by substituting the definitions of the a 
and h coefficients in ( |5.7| ) into ( |5.16D ; we find 
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ds 



A4 



4:\tJ \ds \ dT ^ 
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A4 s fdT\ ( TdXi 
-^~T\Ts)y'^^lT 



Ts 
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dT 
ds 



d^K2\ 



3fi;2 



^^4-©(f)^= 



T 



dK2 

It 



ds J 



3s dT 1 
2T~ds^ 2 



(5.17) 



This is in perfect agreement with the relations obtained in |T0[ using the second law of 
thermodynamics. 

5.4 The Entropy Current 

The entropy of our system is given by 



S 



_d_ 



(To log Z) 



The partition function of our system is given by 



(5.1^ 



\ogZ = --J d'x ^,[P,{Toe-nR + To'P2{Toe-nhr + P3{Toe-n{day\ (5.19) 

(we are careful to explicitly keep track of the temperature dependence in the partition 
function, see the equation ( p. 11 ) for a definition of the functions P). The total entropy 
as evaluated from this partition function is 

d 



S 



dTr 



(To log Z) 



IJVd m - Pi)R + T^ - 3A)/.,r + (^3 - Ps)id^)'] 



(5.20) 



To second order in the derivative expansion, the most general symmetry allowed 
entropy current is given by ||T0[ 



where 

1 



u^ 



(5.21) 



+ ^3 ( R^" - -g^^Rj u, + [A4(u.V)e + A^R + A^^R^pu'^u^) 

+ (Siw^ + 529^ + Sga^)^^ + Bi [{Wsfu" + 256^^5] 
+ [eV^^s - P^\VkU^){VaB^)\ + B^Qa^ + Bja.a^" 

The terms above with Ai and A2 as coefficients are total derivative and do con- 
tribute to the total entropy. It follows that Ai and A2 are unconstrained by comparison 
with equilibrium (even though these terms do not pointwise vanish in equilibrium). 
Terms with coefficients yl4, S2, -B3, -Be and Bj vanish on the equilibrium solution. 
Consequently these coefficients are also unconstrained by the considerations of this 
section. The entropy current coefficients that can be are constrained by comparison 
with (15.201) are A^, A^, Aq, Bi, B^ and B^ 

As above, there are two sources for the second order correction to the entropy of 
our system. The su^ part in Jg contributes to the total entropy at second order in 



derivative expansion because of the second order corrections 6u^ to the equihbrium 
velocity u^^ and 6T to the equihbrium temperature. More precisely, if the equilibrium 
temperature and velocity of our system to second order is given by 

T = T(o) + ST = Toe-"" + 6T and u^ = mJ^) + 6u^ = e~^l, 0, 0, 0) + Su" 

then clearly 

SM°|2nd order = 6^'^ ( — j (5T + s6u^ 



Using ( ^.14|) and (|5.4| ) we find 



(^) ^^ = I [^^' ^ + ^^ ^° ^' - ^^ (^^)' - 2^^ ^'^] 



[P{ i? + P^ ToV' + Ps {daf - 2V. (P3VV)] (5.22) 



s5u' = -2e"'T^ [P^ V^ar + P2 V,/^^*] = -2T^e-"V, {P2 F) 



Therefore using ^l2| , the second order correction to J^, from the perfect fluid piece sm°, 
evaluates to 

S^°|2nd order = ^ [i^l' R + {P2 ' 2^2) T^/' + i^3 {^af] + C-'^V, [2T2p2 /^'a, - P3VV] 

(5.23) 

The second source of two derivative corrections to the entropy current come from 
the explicit two derivative corrections to the entropy current ( ^.2ip evaluated on the 
perfect fluid equilibrium configurations. Using 

o 

R = R- 2{daf - 2VV + —f 

Raf^u'^u^ = {daf + VV + ^f (5.24) 

^o = ^[V,F + 3(V,a)F] 



i?Q — — ( e '^-Roo + flj-Ro 



we find that the zero component of Jg evaluates on equilibrium to 



J°=e- 



A,\Rl-^]+ A,R + A,{Rooe-'n + B,u' + B,{dsf + e" (^ ^ ' '"^ 



uj''\diT) 
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dry dT^ ^ 



^a n2 , ^ ( dB^ A3 rfAs . 



2 V rfT T dT 



-Vj {A,e^''a,r) + V, [{A, - 2A5) VV] 



(5.25) 



Summing ( |5.23| ) and ( p.25| ) and ignoring total derivatives, we find our final result 
for the two derivative correction to the total entropy. 



Total Entropy 



d^xJgi 



A, 



22' 



2A5 + 2^6 - ^3 - 2Si + T2(4P^ - SP^)e-^- 



2a r2 



e^V 
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dT dT^ ^'2 



{day + - 
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dT T 



dM 
dT 



(5.26) 



While the first three terms in ( p.26| ) are Kaluza Klein gauge invariant, the last term 
is not. Let us pause, for a moment to explain this. In subsubsection |2.3.2| we have 
demonstrated that the integral / ^y—giJg is Kaluza Klein gauge invariant provided 
that d^Jg = 0. Now it must certainly be true that the correct entropy current is 
divergence free in equilibrium. However the most general entropy current ( |5.21 



IS 



not divergence free in equilibrium. The non gauge invariant term ion ( ^.26| ) results 
from such terms. The coefficients of these terms must immediately be set to zero 
(even without comparison with a particular form of the entropy). The coefficients 
of the remaining three terms in (|5.26|) must be equated with the coefficients of the 
corresponding terms in ( |5.20| ). In net we have four equations which allow us to solve 
for four of the entropy current coefficients, B^, A3, Bi and B4 in terms of the other 



two {Ar, and Aq) and Pi (the coefficients that appear in the partition function ie.the Pi 
)■ 
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5.4.1 Entropy current with non-negative divergence 

Above we have discussed the constraints on the entropy current from comparison with 
the total entropy of our system. In this subsubsection we will discuss the relationship 
between these constraints and those obtained by imposing the requirement of positivity 
of the entropy current. 

In the study of the positivity of the divergence of the entropy current, it turns 
out that some coefficients in the entropy current are determined in terms of transport 
coefficients, while others are left free (more precisely these coefficients are constrained 
by inequalities involving transport coefficients). The determined coefficients turn out 
to be precisely those that multiply terms that are nonvanishing in equilibrium, namely 
^43, A5, Aq, Bi, B4 and B^. The six equations that determine these six parameters 
are 
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T dT 



1 

4 



-A3 + T 



Ba = — 



dni 

It 

dni 



Ki 



K + 2T-^ + T 



,d^Ki 

It^ 



(5.2^ 



The results ( p. 281) satisfy the constraints ( ^.27|) . In order to verify this one plugs in 



explicit results 
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(5.29) 



for the transport coefficients in terms of action parameters into ( p.28| ) and checks that 
the results are consistent with ( ^.27 ) 



Our results ( |5.27| ) are compatible with but weaker than ( ^■28p . ( |5.28| ) is equiv- 
alent to ( ^.27|) together with A^^ = Aq = 0. As A^ and A^ multiply terms that are 
nonvanishing in equilibrium, we find it surprising that ' our equilibrium study has not 
been powerful enough to demonstrate that A^ and A^ must actually vanish. It is pos- 
sible that we have overlooked a simple principle that forces these coefficients to vanish 
without invoking the principle of entropy increase. 

5.5 The conformal limit 

^^Let us consider Weyl transformation of the full four dimensional metric 

In this subsection ffist we would like to write an partition function which is invariant 
under this transformation. In order to have conformal invariance this partition function 
will have fewer coefficients than the partition function given in ( ^.11|) . Then we shall 



analyze how it will constrain the stress tensor for a conformal fluid. 

Under this transformation several three dimensional quantities transform as fol- 
lows. 

a = a + (p, ai = at, gij = e^'^gij 
(Va)2 = e-2<^ [(Va)2 + 2(Va).(V0) + (V^)^] 

R = e-^^[R-AV^(f)-2{V(f)f] (5.30) 



^z = e'^^^ 
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This subsection has been worked out in collaboration with R. Loganayagam. 



Using ( |5.30| ) we can see that under this transformation the partition function (as 
given in (|5.11|) ) will be invariant (assuming that the total derivative terms will integrate 
to zero) only if the coefficients Pj's satisfy the following constraints. 



Pi(a) = elToe-^ P^{a) = -^ and P^{a) = 2Pi{a) (5.31) 



where Ci and 62 are two dimensionless constants. 
Substituting i^^ in ( |0^ ) we find 



6 = ^4 = C2 = Cs = A4 = 

^2 = 2ki = 2eiToe"'" (5.32) 

A3 = 4X06-^^(61 -262) 

These relations precisely match with our expectation for the independent transport 
coefficients of a conformally covariant stress tensor. Since for a conformally covariant 
stress tensor only two terms {u(^^a^"'u) with coefficient A3 and [R(^u) + -^{^i^}] with 
coefficient ki) can be non zero in equilibrium and a conformally invariant action also 
has only two free parameters, it follows that the existence of a partition function does 
not constrain the stress tensor of a conformal fluid. 

6. Counting for second order charged fluids in 3+1 dimensions 

In this section we will use the methods developed in this paper to answer the following 
question: how many transport coefficients are needed to specify the fluid dynamics of a 
relativistic charged fluid that may not preserve parity, at second order in the derivative 
expansion? We do not attempt to derive the detailed form of the equations so obtained; 
our presentation is merely at the level of counting. If the conjecture at the heart of this 
paper is correct, then an analysis of entropy positivity would yield the same number 
of transport coefficients; however that analysis is much more difficult to perform (even 
at the level of counting), and we do not attempt it here. 

6.1 Parity Invariant case 

Let us first consider the parity invariant case. Table |6.1| list the number of all the the 
parity preserving fluid plus background onshell independent data at second order. From 
this table this we see that the total number of symmetry allowed transport coefficients 
in stress-energy tensor and charge current in landau frame is 

tensors(16) + scalars(18) + vectors(17) = 51. (6.1) 



Type 


fluid+background data 


In equilibrium 


scalars 


16 


9 


vectors 


17 


6 


tensors 


18 


9 



Table 9: parity even data for charged fluids at second order 



Type 


fluid+background 


In equilibrium 


pseudo scalars 


6 


4 


pseudo vectors 


9 


2 


pseudo tensors 


12 


6 



Table 10: parity odd data for charged fluid at second order 



Now let us consider the equilibrium of this system. The third column of table O also 
list the number of scalars, vectors and tensors that can be constructed out of a, Aq, at, 
Ai and g'^^ . The coefficient of these terms that are survive in equilibrium we refer to as 
'non dissipative' coefficients while the remaining we refer to as 'dissipative' coefficients. 
In this case we have a total of 24 non dissipative coefficients. Now there are 9 scalars 
than can be constructed in equilibrium. We list them below 



Hi, V'aV.a, hP, 



p. . pij p. . fij 



VVV.A, VMoV.Ao, V^V.a, V^V.Ao (6.2) 



The last two scalars are total derivatives and hence do not appear in the partition 
function. This tell us that the 24 non dissipative coefficients are determined in term 
of 7 independent coefficients that appear in the partition function which means that 
there will be 17 relation among the 24 non dissipative coefficients. 

In summary the methods developed in this paper predict that parity invariant 
charged fluid dynamics is characterized by 7 non dissipative transport coefficients, 
together with 28 dissipative coefficients (7 scalars, 12 vectors and 9 tensors). Each of 
these 35 transport coefficients is an unspecified function of T and fi. 



6.2 Parity Violating case 

Let us now consider the parity non invariant charged fluids at second order. Table |6.2| 
lists all the parity odd data at second order. From this table we see that number of 
transport coefficients in the parity odd sector is 



pseudo tensors(12) + pseudo scalars(6) + pseudo vectors(9) = 27. (6.3) 



The third column of table |6.2| that out of the 28 parity odd transport coefficients 12 
are non dissipative. 



Now we can construct four new scalars (pseudo) out of the sources. These are 
hsted below 

e'^'d.af.k, e*^'a,Ao/,fc , e'^%aF,k , e''^d,A,F,k (6.4) 

As such all of these scalars listed in ( |6.4D are total derivatives by themselves but only 
two of them can actually be written as total derivatives in the partition function since 
the coefficients that they will appear with are arbitrary functions of a and ^40. 

Any linear combination of the first two scalars can be rearranged as a total deriva- 
tive term and another term that can not be written as total derivative. 

where 



K= f daKi{Ao,a) 



Similar manipulation can be done for the last two scalars listed in 

Thus we see that the 12 parity odd non dissipative coefficients are determined 

in terms of 2 parity odd coefficients in the partition function which means that their 

would be 10 relations in parity odd sector. 

In summary we predict that, at second order, we have 4 parity odd nondissipative 

transport coefficients, together with 2 pseudo scalar, 7 pseudo vector and 6 pseudo 

tensor dissipative coefficients, and total of 20 new coefficients. 

7. Open Questions 

The results reported in this paper suggest several natural follow up questions. In this 
section we list and discuss some of these questions, leaving an attempt to answer them 
to future work. 

The main result of our paper is that two apparently different physical requirements, 
namely the requirement of existence of equilibrium in appropriate circumstances and 
the requirement of the existence of a point wise positive divergence entropy current, give 
the same constraints ^^ on the equations of hydrodynamics in three specific contexts. 
Two questions immediately suggest themselves. Do the results of our paper extend to 
arbitrary order in the derivative expansion, as we have conjectured in this paper? If so. 
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We ignore the inequalities that follow from the principle of entropy increase in this statement. 



why is this the case? Definitive answers to these questions would be very interesting. A 
proof that the existence of equihbrium plus certain inequalities imply the existence of a 
positive divergence entropy current could demystify arguments based on the existence 
of an entropy current, and lead towards a fuller understanding of the second law of 
thermodynamics. 

In the main text of this paper we have derived constraints on the constitutive re- 
lations of hydrodynamics starting from the assumption of the existence of a partition 
function. In the appendices to this paper we have, however, demonstrated that all 
the constraints derived in this paper may also be derived from the weaker assumption 
that fluid admit stationary equilibrium conflgurations in stationary backgrounds. The 
integrability conditions from the demand that the currents and stress tensors in equi- 
librium follow from an action turned out to be automatic in the three examples studied 
in this paper. Is this always the case (we flnd this unlikely). In appropriate situations, 
do the Onsager relations follow from the demand that equilibrium is generated from a 
partition function? 

Hydrodynamical equations that do not obey the constraints described in this pa- 
per are unphysical. Nonetheless these equations are well posed as partial differential 
equations. As a mathematical curiosity one could study such equations in their own 
right. What is the long time behaviour of a generic fluid flow with such equations? 
As equilibrium does not generically exist, flows must, presumably, continue to slosh 
around in a roughly oscillatory manner in the long time limit. It seems clear that such 
indeflnite oscillatory behaviour is inconsistent with the second law of thermodynamics, 
in agreement with the results and conjectures of this paper^^. 

In another direction, the analysis of this paper has led to the consideration of 
partition functions dual to equilibrium hydrodynamics as a function of background 
metrics and gauge flelds. Given a partition function as a function of sources, it is 
standard in quantum fleld theory to Legendre transform this object in order to obtain 
an offshell IPI effective action for the theory. It may be possible implement this 
procedure on our partition function to obtain an offshell action for fluid dynamics 
(albeit one applies only to equilibrium conflgurations) (see [^ for related work). If so, 
what is the interpretation of this action in the context of the fluid gravity map of the 
AdS/CFT correspondence? 

We flnd it interesting that we have been able to encapsulate the effect of the chiral 
anomaly in a 3 + 1 dimensional fluid in a local contribution to the partition function. 
It should be possible to generalize this term to capture the effects of anomalies in 2n 
dimensions for an arbitrary integer n and thereby reproduce the results of [Q. 

It would be very interesting to generalize the work presented in this paper away 
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We thank T. Takayanagi for discussions on this issue. 



from equilibrium. Time dependent partition functions are not in general local func- 
tionals of their sources. These partition functions are, however usually generated by 
coupling local field theory dynamics to sources. Can time dependent correlators (per- 
haps in a Schwinger - Keldysh set up) be generated by minimally coupling a local 
'action' for hydrodynamics to the background metric or gauge field? How does this tie 
in with the fluid gravity map of the AdS/CFT correspondence? 

Apart from the traditional requirement of positivity of the entropy current, and 
the requirement of the existence of equilibrium, emphasized in this paper, one may also 
attempt to constrain the equations of fluid dynamics by demanding that correlation 
functions computed from these equations obey all the symmetry properties that follow 
from the existence of an underlying action (see e.g. 0]). Any system that posseses a 
well defined partition function, as studied in this paper, automatically obeys all these 
symmetry properties for time independent correlators. Do the constraints on hydrody- 
namics that follow from the existence of an equilibrium partition function automatically 
also guarantee that the symmetry requirements on time dependent correlators are also 
met? 

It would be interesting to investigate whether the methods and results of this 
paper carry over to the study of relativistic superfluid hydrodynamics, and in particular 
whether they can be used to rederive the results of [|, |^, |2^, ^, for the most general 
allowed form of the equations of superfluid hydrodynamics at first order. 

Finally, it would be very interesting to investigate the interplay of the principal 
constraint described in this paper (namely the existence of equilibrium for an arbitrary 
static metric) with the AdS/CFT correspondence. Is this constraint merely from the 
structure of AdS/CFT, for an arbitrary bulk Lagrangian, or does it impose constraints 
on possible a' corrections to the equations of Einstein gravity? Within gravity can 
one prove directly that the existence of equilibrium implies the existence of a Wald 
entropy increase theorem (and so the existence [^ of a positive divergence entropy 
current) (see p6[ for related discussion)? 



We feel that several of these questions touch on very interesting, and sometimes 
deep issues in the study of both dissipative systems as well as gravity. We hope to 
report on progress on some of these issues in the future. 
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A. First order charged fluid dynamics from equilibrium in 3+1 
dimensions 

In this appendix we shall rederive the results obtained in section ^ making fewer as- 
sumptions than in that section. In this Appendix we make no reference to the equilib- 
rium partition function, and nowhere assume its existence. The only demand that we 
make on our system is that it admit an equilibrium solution in an arbitrary background 
of the form ( |1 . 1| ) , ( |1.9| ). We also assume that the zeroth order equilibrium configuration 



is given by equation 1.21 



As discussed in section ^ for the parity violating first order charged fluid, one can 
not construct any scalar or pseudo scalar at equilibrium and hence temperature or 
chemical potential does not get corrected to this order. To the first order, the velocity 
corrections can be written as 

5u' = -^e^^V.fe + b2B'j, + b,d'a + hd'Ao (A.l) 

The dissipative part of the stress tensor and the current are written in equation 
37i\ . Since, daj^ — ^, 0, a^y evaluate to zero on equilibrium, we are left with 



Tt^"" 



(A.2) 



We shall now impose that the equations A.2| , A.l obeys the conservation laws 



where 



^ . (A.3) 

V^J^ = CE.B 



2 ^A.4) 



For computational simplicity, we shall take thermodynamic variables temperature T 
and z/ = ;^ as the independent ones. Some useful formulas that are used in computation 



are 



dP 



dT 



V^cu'^ 



T 

e + P dP_ 

dv 



T 
2 



T 



qT 



(A.5) 



uj^W^T, V^B^ 



ji-^ij- 



-2uj^E^ - -B^WT. 



Now using formulas in equation |A.5| , it is straight forward to evaluate the scalar equa- 
tions namely 



V^> = CE.B. 
On setting coefficients of independent data in |A.6| , one obtains 

Tdriiu. + qhi) = 2(e. + qhi), Tdriis + g&2) = {is + 9^2) 
d,{i^ + qhi) = 2T{^B + qb2), d,{^B + 962) = CT 



(A.6) 



d,[{e + P)b2]=T{CB + qb2), drlie + P)b2 



T 



'e + P)b2 



(A.7) 



d,[{e + F)6i] = T{i^ + g6i) + 2T{e + P)b2, drl^e + P)6i] = -(e + P)6i 

ttl = 0^2 = &3 = &4 = 0. 

The vector equation V^^V pT^'^ = V^^J^'^^Jx gives only one new constraint, which is 
given by 

2*, = K|±|l). (A,8) 

On solving |A.7| , one obtains solution for b's and ^'s but with four arbitrary con- 
stants. Now using [A.8| one can eliminate one of the constants in terms of other one. 
Finally we obtain 



bi 



T^ 2 

{-ly^C + 2u^zo + Avz2 + zi) , 



and 



02 = CuT{l - 



e + P'3 

jp2 I 

b2 = '^r^io^'^^ + ^^0 + Z2) 



-^^vT) + T^ [{2vz^ + 22:2) - -^{2v''zo + Avz2 + ^i 
3(e-|-F) '- e + P 

qT 



(A.9) 



-z/T) +T[zfi r^('^^o + ^2)) • 



2(e + P)"'^ ' ' ^'" e + P 

Now identifying zq = 2Co, 22 = C2 and zi = 4Ci we see that equations [A.9| , [A. 10 
exactly same as equations p.ll|, p. 12. 



(A.IO) 



are 



B. First order parity odd charged fluid dynamics from equilib- 
rium in 2+1 dimension 

In this section we shall derive the constraints on parity odd charged fluid dynamics in 
2+1 dimension at first order using just the assumption that there exists a equilibrium 
solution. As discussed in §^ in this case there are 4 transport coefficients and there 
are 6 corrections to the fluid fields. In §^ we were able to express all these 10 functions 
in terms of 2 arbitrary functions in the action ( |4.7| ). This implies that among these 10 
functions only 2 are independent which in turn implies there should exist 8 relations 
among these 10 functions. In this appendix we shall present these 8 relations which 
follows just by demanding that there exists a equilibrium solution. 

We consider the corrections to the fluid fields as in (|4.4| ) and write down the first 
order corrections to the stress tensor and the charge current. The equation of motion 
of fluid dynamics are given by 

(B.l) 
V,J^ = 

Note in particular that the charge current is conserved even in the presence of a back- 
ground gauge field due to the absence of any anomaly in 2+1 dimension. 

Now the scalar equations V^J'^ = and Uy^ (V^T^^ - J^^^Jx) = yields the 
following constraints respectively 

Tj— {xe + P^e) - T—- {xt + pir) + Ptt- {Xe + P^e) = 0, 

a a a ^^.2) 

T— ((e + P)^e) - T— ((e + P)Ct) + ^^ ((e + P)Ce) + T{xt + pir) = 0. 

The vector fluid equations Ppu (V^T^'' — J^^^Jx) = 0, yields the rest of the 6 con- 
straints 

ap dP 

XB = -^rs + ^m„ 

dP dP 

Xn = ^rn + ^mo, 

/ Ti\ /- 9p ap 
, ap ap 

T [XT + Pix) = -^TB + -^niB. 



Note that the first two constraints in ( |B.3| ) are identical to the constraints ( ^.9| ) obtained 



from comparison with most general equilibrium action in § [4 .31 . It is straightforward to 
show that the rest of the constraints in ( [B.2| ) and ( [B.3| ) are solved by the ( |4.1CI| ), ( [4. Ill ), 
(|12D and (|413|) . 



C. Second order uncharged fluid dynamics from equilibrium in 
3+1 dimensions 

In this appendix we will do a similar computation as done in last two appendices for 
3+1 dimensional uncharged fluids at second order. The non-trivial second order (stress 
tensor conservation equation orthogonal to fluid velocity) equation is, 

V,{T,en + e'^ (j2v^V^ + vv) h, + ;^V^nf = (C.l) 

Since, temperature correction is a scalar, we can assume the most generalized 
temperature correction to be of the following form, 

4 

Four dimensional divergence can be expressed as, 

V-^4 V / 2 






= {aA - raTaA)V"^anl + aAV^^n^, (C.3) 

where, we have expressed the two derivative correction to equilibrium stress tensor 



n of ^]^ in a compact form as, (IIjj = aA^fj^ A = 1,11). Using following simple 



derivative formulae 



^'{if^'fjk + ^-^9^k)e'-) = e'-{{y'h,)f\ + 2VV(/,7,, + -^g,,)), (C.4) 



we solve for complete equilibrium solution. In the equation |C.1| , we get following three 
different kinds of terms in parity even sector 

V*(Tensor)jfc, Vfc(Scalar), VfcO"(Scalar), (C.5) 



and following four kinds of terms in the parity odd sector, 



rnkl r _c rj rmn\y2 ^2 rmn 



e™"'V,(V„a/„0, e™"'V,aV„a/„; (C.6) 

Setting the coefficients of V^ (Scalar) to zero, we get the temperature correction as ^^, 
, _ hi 62 , _ &4 , _ h + \{a2-TdTa2) 

rx -fr J^T J^T 



Setting the coefficients of the other terms to zero and using p.7| , we get, the velocity 
corrections as 

(C.8) 



804 — TdxCiA 

Vi = — — , V2 



T^Pt 



04 



and the relations among the transport coefficients as given in |5.16| . Similarly, setting 
the coefficients of the independent terms in the parity odd sector to zero, we get all 
parity odd coefficients zero, that is 



t = V = 61 



S, = 0. 
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